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ABSTRACT 

Theoretical  and  experimental  studies  of  heat  addition  to  external 
supersonic  streams  are  reviewed  following  a  brief  explanation  of  the  basic 
fluid  mechanical  model  and  possible  applications  of  external  burning.  The 
theoretical  section  begins  with  a  terse  review  of  the  extension  of  the  method 
of  characteristics  to  diabatic  flows  and  then  discusses  the  linearized  heat 
addition  models  including  a  new  simplified  method  for  obtaining  the  linearized 
equations.  The  numerous  analyses  of  combustion  via  a  stationary  detonation 
wave  are  categorized  into  four  models  for  planar  heat  addition  and  the  equations 
are  developed  for  the  most  interesting  case  of  the  oblique  Chapman- Jouget 
detonation.  Performance  estimates  are  presented  for  constant  area  and  constant 
pressure  heat  addition  processes  with  and  without  a  consideration  of  the 
expansion  zone  following  heat  addition.  The  governing  equations  for  the  constant 
pressure  analysis  are  developed  for  heat  addition  adjacent  to  both  a  double 
wedge  and  a  flat  plate  surface. 

In  the  experimental  section  all  of  the  available  results  from 
external  burning  tests  are  discussed  beginning  with  the  pioneering  tests  at  NACA 
and  Texaco,  Experiment,  Inc. ,  and  including  tests  of  two-dimensional  and 
axisymmetric  bodies  at  the  Applied  Physics  Laboratory  and  the  Boeing  Co. 

Pertinent  conclusions  ate  drawn  from  each  of  the  tests  and  a  final  compilation 
and  summarization  of  all  the  data  are  given.  It  is  concluded  that  the  maximum 
expected  pressure  coefficient  will  be  near  that  associated  with  a  separated 
boundary  layer  and  that  to  obtain  the  theoretical  maximum  specific  impulse, 
highly  reactive  fuels  with  a  combustion  length  of  a  few  feet  will  be  required. 
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EXTERNAL  BURNING  IN  SUPERSONIC  STREAMS 
Frederick  S.  Billig 
1.  INTRODUCTION 

The  term  "external  burning  ramjet"  (and  the  abbreviation,  ERJ) 
refers  herein  to  any  system  in  which  a  combustible  liquid  or  gas  is 
injected  from  a  vehicle  into  the  external  flow  field  of  the  body  and 
burns,  thereby  altering  the  flow  field  and  producing  forces  on  the  body 
due  to  the  combined  effects  of  flow  interactions  and  heat  addition. 

Although  subsonic  applications  of  external  burning  may  be  feasible,  only 
supersonic  flight  velocities  are  considered.  Systems  of  this  type  have 
been  rather  extensively  studied  during  the  past  several  years  and  experi¬ 
ments  have  been  made  which  have  demonstrated  the  feasibility  of  the  external 
burning  concept,  but  in  general  they  have  shown  only  part  of  the  performance 
potential  of  the  system.  The  available  unclassified  literature  in  this  area 
is  reviewed  and  conclusions  are  drawn  regarding  the  limitations  of  external 
burning  systems  based  principally  on  the  phenomena  observed  in  the  tests 
reported. 

The  generation  of  useful  forces  on. the  external  surfaces  of  an 
aerodynamic  body  requires  deflection  of  the  streamlines  in  the  flow  field 
about  the  body  in  such  a  manner  that  increased  pressures  are  produced  on 
chosen  surface  areas.  Consider  the  external  flow  fields  above  the  flat 
plates  sketched  in  Fig,  1.  Undisturbed  streamlines  are  shown  in  (a),  and 
the  other  sketches  show  the  disturbed  streamlines  pattern  caused  by  an 
aerodynamic  flap  (b),  mass  addition,  (c  ) ,  and  mass  addition  plus  heat  addition 
(external  burning)  (d).  The  correr. ponding  surface  pressure  profiles  also  are 
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Fig,  1  STREAMLINE  PATTERNS  FOR  VARIOUS  FLAT  PLATE  SYSTEMS 
IN  A  SUPERSONIC  FLOW,  AND  THE  CORRESPONDING  SURFACE 
PRESSURE  PROFILES 
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shown  in  Fig.  1  by  che  correspondingly  lettered  curves.  For  case  (b),  the  pres¬ 
sure  rise  due  to  streamline  deflection  at  station  1  must  be  followed  by  expan¬ 
sion  to  below  ambient  pressure  at  2  and  then  recompression  to  about  p^  at  3.  Thus, 
to  obtain  the  greatest  net  positivi  force,  the  flap  should  be  positioned  suffi¬ 
ciently  far  aft  so  that  station  2  corresponds  to  the  trailing  edge  of  the  plate. 
Furthermore,  there  is  an  attendant  drag  force  penalty  on  the  flap.  For  sim¬ 
plicity,  in  cases  c  and  d  the  heat  and/or  mass  addition  is  confined  to  stream- 
tube  I  in  a  zone  of  finite  length  1-2  and  is  assumed  to  occur  at  constant 
pressure.  The  surface  pressure  is  sustained,  however,  until  the  first  expansion 
wave  strikes  the  surface  at  3.  then  the  pressure  declines  to  a  value  near  p^  at  4. 
In  effect,  the  heat  addition  (and  the  mass  addition  to  a  lesser  degree)  represents 
a  volume  source  in  streamtube  I  which  turns  adjacent  streamtubes  (II,  III,  etc.) 
gi'  ng  a  pressure  rise,  similar  to  the  case  of  the  flap  but  *?ith  a  significantly 
lesser  expansion  effect  and  no  drag  penalty.  In  addition  to  the  pressure  force 
there  is  a  reaction  force  caused  by  injection,  which  has  components  in  the 
thrust  and/or  lateral  directions,  depending  on  the  angle  of  injection. 

The  foregoing  crude  description  of  the  effects  of  the  heat  and  mass 
addition  zones  is  oversimplified,  because  they  need  not  be  zones  of  constant 
pressure,  and  the  details  of  boundary  layers  and  possible  attendant  separated 
zones  have  been  omitted.  However,  from  sketch  d  it  is  apparent  that  to  obtain 
the  greatest  total  normal  force  from  the  positive  pressure  field  developed  b> 
the  external  burning  case,  it  is  necessary  to  extend  the  surface  to  the  end 
of  the  expansion  zone  (point  4).  On  the  other  hand,  a  higher  force  coefficient 
f force/ (dynamic  pressure  x  area)]  would  be  obtained  if  the  plate  were  cut  off 
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at  point  3.  This  point  and  the  effect  of  initial  pressure  level  are  discussed 
further  in  a  later  section. 

The  possible  applications  of  the  external  burning  principle  fall 
into  three  classes:  (A)  side-force  generating  devices  for  attitude  control, 
(B)  thrust-generating  (or  drag-reducing)  devices,  and  (C)  devices  which 
produce  both  thrust  and  attitude  control  (or  lift).  Sketches  of  possible 
configurations  are  shown  in  Fig.  2.  The  attitude  controller  for  an  axi- 
symmetric  vehicle  [Fig.  2(A)]  has  injection  aft  of  the  center-of-gr avity 
in  any  one  of  four  quadrants.  Longitudinal  "fences"  separate  the  quadrants 
to  reduce  the  dissipation  of  the  positive  pressure  field  through  circum¬ 
ferential  spillover.  Note  that  the  downward  force  due  to  external  burning 
leads  to  positive  pitch  and  therefore  puts  the  external  burning  region  in 
the  leeward  zone,  which  could,  at  large  pitch  angles,  produce  adverse  con¬ 
ditions  for  combustion.  However,  if  the  external  burning  is  being  used 
solely  to  trim  the  body,  then  it  is  conceivable  to  design  an  aerodynami” 
cally  unsteady  vehicle,  in  which  case  the  external  burning  will  always  occur 
in  the  windward  zone.  Attitude  control  systems  based  on  external  burning 
ahead  of  the  c.g.  are  conceivable  but  appear  to  be  less  attractive  due 
to  the  difficulty  of  confining  the  positive  pressure  field  to  produce 
an  effective  pitching  moment.  The  thrust  generating  device,  Fig.  2(B),  could 
conceivably  be  either  the  total  vehicle  or  a  podded  or  airfoil  engine.  At 
the  "knee"  fuel  is  added  to  the  air  compressed  by  oblique  shock  and/or 
isentropic  turning  on  the  forebody,  and  combustion  maintains  a  positive 
pressure  field  on  the  aft  body  which  is  greater  than  that  or  the  compression 
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surface,  thus  producing  net  thrust.  Both  axisymmetric  and  two-dimensional 
configurations  are  possible  for  this  case  and  the  combined  thrust  and  lift 
case,  2(C).  The  combined  device,  which  has  a  flat  top  (hence  no  positive 
pressure  on  top)  but  a  positive  pressure  field  over  its  entire  lower  surface 

develops  considerable  lift.  It  could  be  used  as  a  "propulsive  wing"  or 

"external  burning  ramjet  (ERJ)".  Analysis  of  configurations  (B)  or  (C) 
has  shown  them  to  have  efficient  thermodynamic  cycles  only  at  very  low 
thrust  levels,  i.e.,in  cruising  flight,  as  shown  later.  For  accelerating 
missions,  ducted  conventional  ramjets  or  supersonic  combustion  ramjets 
(scramjets)  are  considerably  more  efficient. 

2.  LITERATURE  REVIEW  OF  THEORETICAL  ANALYSES 

_ OF  EXTERNAL  BURNING _ 

Theoretical  analyses  of  the  effect  of  heat  addition  to  external 
supersonic  flow  fields  began  to  appear  in  the  literature  about  fifteen  years 
ago,  and  the  first  (then)  classified  experimental  results  followed  one  to  two 
years  later.  Following  the  very  fundamental  work  in  one-dimensional  diabatic 
flow,  ^  ^  Pinkel  and  Serafini^  extended  the  method  of  characteristics  to 
include  the  effects  of  heat  addition  in  an  irrotational  supersonic  flow  and 
developed  a  graphical  method  of  solution  for  shock-free  flow  with  heat  addition. 
Using  this  technique  it  is  possible  to  find  an  exact  solution  (with  the  above 
constraints  to  the  flow  field  and  pressure  distribution)  for  flows  having 
continuous  total  temperature  variation  in  the  streamwise  direction.  In 
Reference  11  this  method  was  used  to  determine  the  pressure  distribution  and 
aerodynamic  coefficients  of  a  symmetrical  circular  arc/wing  (Fig.  3)  for  flight 
Mach  numbers  (Mq)  of  3  and  5.  For  moderate  total  temperature  ratios  over 
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the  heat  addition  zone  (T  /T  =  1.243  and  1.126),  the  heat  release  beneath 

t2  tL 

the  middle  portion  of  the  airfoil  produces  a  significant  pressure  rise, 
markedly  increasing  the  lift  and  slightly  decreasing  net  drag  (5-107.)  so  that 
the  lift/drag  ratio  (L/D)  and  quarter-chord  moment  coefficient  are  increased 
by  factors  of  1.7  to  2.2,  and  2.5  to  4.2,  respectively,  compared  to  no-heat- 
addition  values. 

One  of  two  motivations  has  inspired  succeeding  authors  to  find  other 
methods  of  analyzing  external  burning: 

1)  Greater  simplicity. --The  method  of  Ref.  10  is  tedious  and 
does  not  lend  itself  to  simple  evaluation  of  the  important  variables  which 
affect  performance.  The  approaches  used  to  simplify  the  problem  involve  either 
linearization  of  the  equations  of  motion  or  postulation  of  a  one-dimensional 
heat  addition  process  (e.g.,  constant  area  or  constant  pressure). 

2)  Greater  realism. — As  experimental  data  from  external  burning 
tests  have  accumulated  it  has  become  possible  to  postulate  new  analytical 
models  which  include  effects  that  have  been  observed  and  lead  to  closer 
correlation  of  theory  with  experiment. 

2.1  Linearized  Solutions 

An  approximate  formula  for  a  linearized  solution  for  the  pressure  field 

generated  by  a  moderate  rate  of  heat  release  was  developed  from  first  principles 
12  13 

by  Chu.  Gazley  arrived  at  the  identical  result  by  employing  a  "piston" 
concept  to  translate  the  effect  of  heat  addition  into  an  effective  deflection 
of  the  flow.  The  result  given  by  these  authors  can  also  be  obtained  by 
considering  the  streamline  deflection  caused  by  heat  addition  in  the  simple 
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model  shcvn  in  Fig.  4(A).  Heat  is  added  co  streamtube  I  in  the  region 

from  1  to  2,  increasing  the  streamtube  height  from  Yj  to  Y9  and  causing 

a  deflection,  5„,  to  adjacent  streamlines.  If  a  is  the  rate  of  heat  addi- 
n 

tion  per  unit  area  and  Q  is  the  heat  addition  per  pound  of  air,  then  from 
continuity  and  the  equation  of  state 

P,u  Y  =  P2u?Y2,  or  PjUjYj/Tj  =  P^Yj/T,  (1) 

and  from  the  geometry. 


fY  /cos  (6h/2)]  -  Yt 


X, 

t 


=  tan  o 


H 


From  linear  theory  (small  deflections) 


u2~ur 


cos  (6jj/2)  «  1, 


and  tan  6  =»  6 
H  H 


(2) 

(3) 


Therefore 

(Y?  -  Y1)/Xf  ^  oH  (4) 

From  conservation  of  energy, 

T2  *  Tx  +  Qgj/  cp  (5) 

From  conservation  of  momentum,  p^  %  p.,  since  «  u^;  Eq.  (1)  becomes 

VT1  =  Y2/T2  (6) 

therefore,  substituting  (5)  and  (6)  into  (4), 


hL 

<2£j' 

-  Y 

Ti  1 

1  CP| 

1  Qg-T  \ 

(7) 

CP 

T1  Xf 

cPTlxf 

QgJ  Yj 

P1  Ul/X 

f 

(8) 
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A)  GENERAL  THEORETICAL  MODEL 


B)  FLAT -PLATE  AIRFOIL  WITH  LINEARIZED  HEAT 
ADDITION  ON  LOWER  SURFACE 


Fig.  4  LINEARIZED  HEAT  ADDITION  MODEL 
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SO 

resulting  in  the  same  expression  for  the  turning  angle,  6^,  as  a  function 
of  the  ratio  of  heat  addition  per  unit  area  as  was  developed  in  Refs.  12 
and  13.  Pressure  coefficients  for  various  heat  additions  can  then  be  obtained 
by  using  the  familiar  linearized  supersonic  flow  result, 

C  =  2  b/ (M  2  -  1)^  (10) 

P  0 

Since  linear  solutions  may  be  superimposed,  a  general  equation  for  the 
pressure  coefficient  of  a  body  with  heat  addition  deflection  6^  and  surface 
inclination  6^,  can  be.  obtained,  viz: 

C  =  2  (5  +  6  )/(M  2  -  1)^  (ID 

p  bHu 

Using  this  approach  Gazley  showed  that  the  lift  and  drag  coefficients  for 
a  flat  plate  airfoil  at  angle-of-attack,  or,  with  heat  addition  over  a  fraction, 
f  =  X^/C,  of  the  plate  (Fig.  4B)  are: 

CL  =  4  (a  +  £  or  f  6h)/(Mq2  -  1)?  (12) 

and 

CD  =  4  {a2  +  £  a  f  6h)/(Mq2  -  1)^  (13) 

For  a  biconvex  airfoil  of  thickness-chord  ratio  of  t/C  the  same  expression 
Eq.  (11)  holds  for  and 


CP  °iuiTi 


%  piui/R 


Y  PjU, 
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He  then  shows  (Fig.  5)  for  the  5-percent- thick  biconvex  airfoil  that 
the  linear  and  exact  solutions  of  Pinkel,^  et  al  agree  well  at  low  heat 
addition  rates  and  start  to  deviate  at  higher  heating  rates  where  both 
methods  become  inaccurate.  The  calculations  show  that  for  small  negative 
angles-of-attack  and  heating  deflections  of  greater  than  5  or  6  degrees, 
net  thrust  can  be  produced.  Moreover,  by  defining  a  specific  impulse  based 
on  (drag  cold  -  drag  hot)/(fuel  flow  rate)  for  a  given  lift  force,  where  the 
angle-of-attack  for  the  hot  flow, 

“h  -  “c  -  i  £  V  (15) 

is  smaller  than  the  cold  value,  specific  impulses  of  300  to  1800  sec  are 

computed  for  a  fuel  with  a  heating  value  of  20,000  Btu/lb. 

14 

Mager  arrived  at  essentially  the  same  results  as  Gazley  and  showed 

good  agreement  with  the  more  exact  solution  of  Ref.  10  for  the  circular  arc- 

airfoil.  He  started  with  the  general  definition  of  the  pressure,  velocity 

and  density  disturbances  caused  by  a  heat  source  in  supersonic  flow  as  sug- 
15 

gested  by  Tsien  and  used  a  slightly  different  method  of  linearization.  In 
addition,  he  developed  a  set  of  performance  indices  to  compare  external 
burning  (ERJ)  configurations  with  a  winged  vehicle  using  conventional  ducted 
subsonic  combustion  ramjecs  (CRJ).  He  concluded  that  the  energy  requirements 
during  cruise  were  comparable  for  the  two  engine  systems.  He  included  the 
estimated  skin  friction  drag  and  concluded  that  it  is  more  profitable  to 
take  advantage  of  the  additional  force  generated  by  heat  addition  by  decreasing 
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Fig.  5  LIFT  OF  SUPERSONIC  AIRFOIL  WITH  HEAT  ADDITION 
ON  LOWER  SURFACE.  (Ref.  13) 
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wing  area  (and  friction  drag)  and  operating  at  the  same  angle-of-attack  rather 
than  maintaining  the  same  wing  area  and  reducing  or  (and  wave  drag) . 

2.2  Planar  Heat  Addition 

Another  simplified  method  of  approach  to  the  theoretical  solution 
of  the  external  heat  addition  problem  is  to  postulate  planar  heat  addition 
at  some  angle  to  the  air  flow,  i.e.,  from  a  physical  standpoint,  an  infinitely 
fast  heat  release.  The  theoretical  models  which  have  been  postulated  include 
both  normal  and  oblique  planar  heat  addition,  with  or  without  accompanying 
shock  waves,  as  shown  in  Fig.  6.  Figure  6(A)  shows  a  simple,  normal,  planar 
heater,  which  in  effect  is  a  one-dimensional,  constant-area  heat  addition. 

All  of  the  changes  in  properties  across  the  heater  can  be  obtained  from  the 
well-known  one-dimensional  equations. ^  For  incoming  supersonic  flow 
(Mq  >  1)  of  streamtube  I  there  is  a  pressure  rise  due  to  head  addition  and 
a  corresponding  downstream  Mach  number,  >  1.  Adjacent  streamtubes 
(II,  III,  etc.)  must  compensate  for  this  pressure  rise  by  turning  through 
an  oblique  shock  wave.  The  incompatability  in  flow  direction  and  pressure 
downstream  is  then  resolved  by  a  series  of  expansion  waves  downstream. 

In  Fig.  6(B)  the  normal  heater  is  preceded  by  a  normal  shock  wave,  so  that 
the  flow  upstream  of  the  heater  (1)  is  subsonic.  Subsonic  constant-area 
combustion  must  then  be  postulated,  with  M2  <  1.  This  process  is  generally 
referred  to  as  a  strong  detonation,  and  its  limiting  case  for  =  1,  when  the 
distance  between  the  shock  wave  and  the  heating  wave  becomes  vanishingly  small, 
is  called  a  Chapman  Jouguet  detonation.  Although  solutions  of  the  equations 
of  motion  permit  both  the  reverse  order  of  the  processes  and  supersonic  Mach 
numbers  at  2,  the  very  unrealistic  kinetic  situation  of  combustion  occurring 
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Fig.  6  THEORETICAL  MODELS  FOR  PLANAR  HEAT  ADDITION 
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at  lower  rather  than  higher  pressure  and  temperature  forbids  the  former, 
and  entropy  considerations  rule  out  the  latter.  Again,  pressure  compensa¬ 
tion  in  the  external  stream  can  be  accomplished  by  an  oblique  wave,  because 
although  the  pressure  at  is  the  normal-shock  value,  the  pressure  at  2 
is  lower  (subsonic,  constant-area  heat  addition),  so  that  a  shock  weaker 
than  a  normal  shock  will  balance  the  pressure.  Rather  than  considering 
the  difficult  expansion  case  of  <  1.0,  most  authors  have  chosen  to  study 
the  Chapman  Jouguet  case  of  =  1.0,  which  for  the  model  shown,  will  only 
occur  at  a  specified  temperature  ratio,  T^/T^ ,  for  each  free  stream  condition. 
To  avoid  this  limitation,  a  subsonic  compression  (streamtube  enlargement) 
of  the  flow  is  postulated  in  region  1'  in  such  a  manner  that  the  condition 
of  =  1.0  is  always  met. 

Figure  6(C)  shows  an  oblique  heater  of  normal  height  Y^.  This  is 
analytically  handled  by  applying  the  one-dimensional,  constant-area  heat 
addition  equations  to  the  normal  component  of  the  upstream  velocity.  Note, 
however,  that  the  downstream  velocity  vector,  u^>  must  be  directed  away  from 
the  surface,  which  means  that  some  type  of  separated  zone  would  have  to  exist 
downstream  of  the  heat  addition  plane  and  would  have  to  extend  to  a  height, 
h2 : 

h2  =  1  -  (P1u1/p2u2)  cos  Pg  (16) 

h2  Sin  (6H  "  V  (17) 

and  the  tangent  of  the  separation  angle  is  tan  6  =  — — - - -  v  1 

S  Yi  "  h2 

Downstream  of  X2  the  flow  would  have  to  expand  and  turn  to  adjust  pressure 
and  flow  direction.  The  flow  in  streamtube  II  which  does  not  undergo  a  rise 
in  stagnation  temperature  needs  to  be  turned  through  a  lesser  angle  (by  an 
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oblique  shock)  Co  reach  than  does  the  flow  in  streamtube  I,  so  that  a 
vortex  sheet  must  exist  downstream  of  X Note  that  in  these  simplified 
models  the  heat  addition  is  limited  to  streamtube  I  and  there  is  no  mixing 
with  adjacent  streamtubes.  Therefore,  even  though  the  pressure  is  matched 
along  the  boundary  between  I  and  II,  the  entropy  is  greater  in  I,  and  a 
vortex  sheet  must  divide  these  two  regions  of  the  flow. 

The  last  case,  shown  in  Fig.  6(D)  and  detailed  in  Fig.  7,  is  the 
strong  oblique  detonation,  i.e.,  the  oblique  planar  heat  addition  is  pre¬ 
ceded  by  an  oblique  shock  wave,  and  the  component  of  velocity  normal 

to  the  heater  plane  angle,  6  ,  is  subsonic.  Again,  the  special  case  of 

=  1.0  would  be  called  an  oblique  Chapman  Jouquet  detonation.  Presumably 

all  combinations  of  oblique  shock  angles,  9,  and  heater  angles,  6  ,  are 

n 

possible  as  long  as;  a)  the  flow  behind  the  oblique  shock  is  supersonic, 

b)  9  is  below  the  shock  detachment  angle  for  the  local  Mach  number,  and 

c)  the  required  heat  release  is  not  excessive.  However,  as  noted  by 

Willmarth, ^  for  each  9  there  is  one  value  of  6  which  results  in  a  velocity 

h 

vector  u^  after  heat  addition  which  is  aligned  to  the  surface.  The  model 
based  on  this  particular  situation  has  the  virtue  of  not  requiring  a  separated 
zone  as  in  Fig.  6(C).  The  vector  diagram  shown  in  Fig.  7  for  this  special 
case  shows  that  for  a  given  9  and  6  the  tangential  components  of  velocity 
are  related 


V  =  U1T  =  U0  C0S  9 


(18) 


(19) 
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Using  the  example  chosen  in  Ref.  17  of  =  2.0  at  30,000  ft  altitude  hut 
extending  the  planar  heat  addition  analysis  to  large  heater  angles  gives 
the  results  shown  in  Figs.  8  and  9.  In  Fig.  8,  q,  the  heating  rate  per  unit 
area  is  shown  as  a  function  of  the  surface  pressure  rise,  p?  -  p^,  for  various 
heater  angles,  5  .  Lines  of  constant  turning,  9  -  6  ,  which  imply  constant 
shock  angles,  9,  are  also  shown.  The  curves  are  bound  on  the  left  by  the 
Chapman  Jouguet  limit  (M^  =  1.0)  at  low  pressure  rise  and  by  the  maximum 
attached-shock  strength  at  higher  pressure.  On  the  right,  the  curves  are 
bounded  by  a  maximum  heating  rate,  because  the  temperature  rise  across  the 
heater  becomes  infinite  as  6  "*  (9  -  o  ).  For  a  given  pressure  rise  the 
oblique  detonation  processes  requi. e  more  heat  release  than  that  predicted 
by  linear  theory  because  of  the  inclusion  of  the  turning  losses.  In  Fig.  9 
this  effect  is  shown  in  curves  of  pressure  rise  per  unit  energy  release,  or 
"lifting  efficiency"  for  various  turning  angles.  As  the  heater  and  turning 
angles  approach  zero,  the  linear  theory  value  of  Ap/q  obtained  by  combining 
Eq.s  (9)  and  (10), 


AR  = 
q 


(28) 


or  4.65  x  10  sec/ft  in  the  example,  is  approached. 

Woolard^’^  studied  the  same  case  of  strong  oblique  detonations  and  lists 

tables  of  properties  for  a  range  of  conditions  in  Ref.  18.  Similar  tables  have 

20 

been  generated  by  Chinitz,  et  al.  In  Ref.  19,  Woolard  develops  a  set  of  equa¬ 
tions  similar  to  Eqs.  (18-27)  for  conical  rather  than  two-dimensional  flows. 
Parametric  curves  are  presented  for  Chapman  Jouguet  detonative  flows  about  cones. 
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Refinements  to  the  oblique  planar  heater  model  of  Wilmarth  can  be  made.  One 
of  these  is  to  consider  multiple  oblique  heater  planes  having  a  common  origin 
rather  than  a  single  planar  heater.  For  example,  in  the  case  of  two  oblique 
heater  planes,  the  analysis  for  the  first  heater  plane  is  the  same  as  before 
except  that  the  flow  velocity  vector  behind  this  heater  is  directed  away  from 
the  plate.  Turning  through  the  second  heater,  which  can  be  described  by  Eqs. 
(18-27),  changing  subscripts,  then  aligns  the  flow  with  the  plate.  Typical 
calculations  made  at  =  6.0  and  Z  =  20,000  ft  showed  that  for  a  given  heater 
flux  the  two  heater  plane  case  predicts  a  3-3.5  per  cent  lower  pressure  rise 
than  the  single  heater.  In  effect,  in  the  limit,  an  infinite  number  of  weak 
centered  planar  heaters  becomes  a  wedge-shaped  continuous  heat  addition  pro¬ 
cess.  From  a  physical  point  of  view,  this  implies  the  rather  untenable  sit¬ 
uation  of  a  reaction  requiring  a  rate  of  temperature  rise  in  the  streamwise 
direction  that  is  proportional  to  the  distance  from  the  surface. 

A  second  extension  is  to  consider  the  expansion  region  downstream  of 
the  constant  pressure  zone  for  a  finite  heater  height  (Y^,  Fig. 7  ).  The 
pressure  decays  iri  the  region  behind  the  first  Mach  line  emanating  from  the 
terminal  point  of  the  planar  heater,  thus,  the  surface  pressure  drops  when 
the  Mach  line  reaches  the  surface.  This  decaying  (but  still  positive)  pres¬ 
sure  field  which  occurs  without  an  additional  expenditure  of  energy  therefore 
results  in  a  higher  specific  impulse  than  that  shown  by  Wilmarth.  Several 
methods  of  calculation  of  this  decaying  pressure  field  have  been  suggested 
by  other  authors  and  are  described  in  the  subsequent  discussion. 
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21 

Luidens  and  Flaherty  applied  the  normal  weak  detonation  model  of 
Fig.  6(A)  to  the  ERJ  model  of  Fig.  10(a).  Calculations  were  made  for 
various  combinations  of  wedge  angles  and  at  different  angles-of-attack 
in  order  to  determine  the  optimum  wing  geometry  and  angle-of-attack.  The 
ERJ  performance  was  based  on  a  specific  lift  parameter  L/h  w^,  where  h  is 
the  heating  value  of  the  fuel,  which  can  be  related  to  range  through  the 
conventional  Breguet  range  equation  for  the  case  when  thrust  equals  drag 


and  lift  equals  weight: 
u 


R  =  T- 
wf 


0 


1  -  u 


V  2 


An  “ 


w, 


(29) 


B 


Most  of  the  calculations  were  made  using  the  linearized  theory  and  showed 
that  "...  a  wing  designed  for  maximum  L/h  w^  will  have  the  following  character¬ 
istics:  a  moderate  thickness  ratio,  a  flat  top  surface,  and  the  maximum¬ 
thickness  point  of  the  wing  well  downstream  (i.e.,  e  >  e^).  It  should  be 

operated  at  maximum  angle-of-attack  (i.e.,  1°  to  3°) . "  More  rigorous 

r  ’culations,  using  the  method  of  characteristics  for  the  downstream  expansion 
fTig.  6(A)  ]  were  made  and  showed  efficiencies  of  about  half  those  estimated 
by  linear  theory.  A  few  calculations  were  made  to  show  that  performance  could 
be  improved  by  distributing  heat  sources  in  Region  II  of  Fig.  10(A).  Finally, 
they  compared  ERJ’s  with  conventional  ramjet  (CRJ)  cruise  vehicles  with  wings 
and  obtained  Fig.  10(B).  Curves  A  and  B  are  for  CRJ's,  A  being  quite  optimistic 
and  B  more  conservative.  Curve  C  represents  optimistic  performance  of  an  ERJ 
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with  a  distributed  heat  source,  and  curve  D  is  for  an  ERJ  with  a  single 
planar  flame.  They  conclude  that  at  Mach  numbers  greater  than  8,  the  ERJ 

has  greater  cruise  efficiency  than  the  (subsonic  combustion)  CRJ  vehicle. 

■?2 

Lomax  also  considered  distributed  heat  sources  throughout  region  II  of 
Fig.  10(A)  and  found  sizeable  gains  (>  200’4)  above  linear  theory,  however, 
it  is  indeed  difficult  to  imagine  how  such  a  fuel  distribution  could  be 
obtained  in  practice. 

Smith  and  Davis  at  Experiment  Inc.  were  the  first  to  consider  the  normal 
strong  detonation  case  of  Fig. 6(B).  They  had  reasoned  that  the  stabilization  of 
a  supersonic  flame  would  be  difficult  and,  therefore,  postulated  a  subsonic 
heat  addition  preceded  by  a  normal  shock.  The  particular  case  of  =  1.0 
was  chosen  for  study  and  the  pressure  distribution  behind  the  planar  heater 
was  found  by  balancing  the  pressure  in  streamtube  I  with  that  in  external 
streamtubes  by  taking  successive  oblique  turns.  In  region  1,  between  the 
shock  and  the  heater,  the  pressure  coefficient  was  based  on  the  linear 
theory  value  [Eq.  (10)]  which  now  introduces  a  geometric  factor,  the  ratio 
of  flame  height  to  shock-flame  separation  distance,  ***  ^or  st0ichi0" 

metric  burning  of  hydrogen  adjacent  to  a  flat  plate  they  obtained  side 
force  specific  impulses  of  2140,  3560  and  4860  secs  at  Mach  numbers  of  2,  3 
and  4.  For  lean  limit  operation,  (no  separation  between  shock  and  heater) 
the  corresponding  impulses  were  7200,  8800  and  9600  sec.  They  also  applied 
the  planar  flame  concept  to  a  double  wedge  (equal  fore  and  aft  wedge  angles) 
model  similar  to  Fig.  2(B)  and  obtained  a  specific  impulse  defined  as  (thrust 
with  burning)  -  (cold  flow  drag)  t  (fuel  flow  rate).  .  Table  I,  their  M 
is  the  local  Mach  number  on  the  aft  wedge  ahead  of  the  normal  shock.  This 


-26- 


THE  JOINS  HOPKINS  UNIVERSITY 

APPLJED  PhiYSICS  LABORATORY 

StL>CM  Sr*iNG  MAMVUAHD 


author  calculated  the  corresponding  flight  Mach  numbers,  Mq,  shown  in 
Table  I,  based  on  a  leading-edge  oblique  shock  and  Prandtl-Meyer  turning 
at  the  "knee". 

23 

Monchick  and  Dugger  considered  the  case  of  strong 
detonation  on  the  aft  surface  of  the  two-dimensional  airfoil  shown  in 
Fig.  11,  which  in  effect  is  the  Smith  and  Davis  model  split  on  the 
centerline.  For  the  particular  case  of  wedge  angles  (o')  of  10°  and  stoic! io- 
metric  hydrogen  combustion  with  sonic  burned  gas,  the  flame  height  to  surface 
length  ratio,  h/it  was  found  for  the  condition  of  zero "net  thrust. 


TABLE  I 


Smith  and  Davis  Results  for  Specific  Impulse  on  Double- Wedge  Model 

Hedge  M  „  x  (Thrust  HotHteag  Cold).  sec 

Half-Angle  Of  w^ 

Lean  Limit 

Stoichiometric 

5° 

2.0 

1.87 

630 

185 

3.0 

2.76 

770 

310 

4.0 

3.65 

840 

425 

10° 

2.0 

1.67 

1250 

370 

3.0 

2.56 

1530 

620 

4.0 

3.38 

1670 

850 

Figure  11(B)  shows  that  h II  would  exceed  0.7b  even  for  flight  Mach  numbers 
in  the  range  of  3  to  5.  This  result  was  rather  discouraging  from  a  practical 
design  standpoint,  and  an  endeavor  was  made  to  make  a  more  refined  analysis 
of  the  expansion  process  to  see  if  the  result  was  due  only  to  the  simplifying 
assumptions.  Since  the  expansion  process  in  region  5  begins  at  the  hot-cold 
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B)  MINIMUM  RATIO  OF  PLANE  FLAME  HEIGHT  TO  EXPANSION  RAMP 
LENGTH  REQUIRED  FOR  NET  THRUST 


C)  COMPARISON  OF  CHARACTERISTICS  WITH  CHANNEL  FLOW  SOLUTION 
FOR  PRESSURE  DISTRIBUTION  ON  SURFACE  BEHIND  FLAME 


Fig.  11  TWO  DIMENSIONAL  AIRFOIL  WITH  NORMAL  PLANAR  HEAT  ADDITION.  (Rsf.  23) 
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interface  (see  Fig.  6)  the  wedge  surface  pressure  decay  in  the  actual 
expansion  process  would  first  lag  and  then  lead  the  interface  pressure 
decay  used  in  the  first  study.  Thus,  the  integrated  net  force  on  the 
aft  wedge  could  result  in  a  larger  thrust  in  the  actual  expansion  case. 

For  the  same  a  =  10°  model,  a  method-of-characteristics  solution  for 
Mq  =  5,  with  (i.e.,  no  subsonic  compression  and  below  stoichio¬ 

metric  fuel-air  ratios)  and  =  1.0  gave  the  surface  pressure  decays 
shown  in  Fig.  11(C)  (solid  curve).  Integration  of  the  pressure  force 
showed  that  the  minimum  h/A  for  the  characteristics  solution  was  0.42 
compared  to  1.4  for  the  "channel  flow"  solution.  From  Fig.  11(B), 

(h/i)  .  is  0.75  at  M.  =  5,  thus  the  M„  =  M.  case  requires  a  larger 
flame  height  than  the  stoichiometric  (ER  =  1.0)  case.  Presumably,  if 
the  characteristics  solution  were  applied  to  the  stoichiowetric  case, 
lower,  and  therefore  more  physically  realizable,  flame  heights  would  occur. 

A  further  refinement  in  the  planar  flame  model  was  suggested  by 
Woolard  to  handle  the  subsonic  compression  region  between  the  normal  shock 
and  the  planar  heater.  Instead  of  assuming  the  linearized  flow  solution 
in  region  1  of  Fig.  12,  an  approximate  but  more  realistic  description  of 
the  flow  was  made.  The  shape  of  the  detached  shock  wave  from  A  to  B  was 
assumed  to  be  a  hyperbola  becoming  assympotic  to  the  free  stream  Mach  line 
C-D  at  point  B.  The  shape  of  the  hyperbola,  attachment  distance  and  the 
average  pressure  coefficient  for  region  l  were  based  oj|  the  method  of  Moeckel. 
As  before,  the  second  shock  wave  emanating  from  the  top  of  the  heater  balances 
the  pressure  between  the  burned  and  unheated  stream  tubes.  The  flow  in 


25 
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streamtube  II  is  analogous  to  the  spillover  flow  in  a  ramjet  inlet 

operating  supercritically.  As  the  mass  flow  in  this  streamtube  approaches 

zero,  Xg  goes  to  zero  and  the  flow  situation  is  again  that  of  Fig.  6(B) 

(critical).  Woolard  also  developed  an  approximate  method  for  obtaining 

the  pressure  decay  downstream  of  the  heater  based  on  an  exponential 

pressure  decay  beyond  the  last  simply  reflected  Prandtl-Meyer  characteristic. 

Charts  for  determining  all  of  the  pertinent  flow  field  parameters  and  the 

* 

normal  force  coefficients  as  a  function  of  the  rate  of  heat  addition  are 

presented  for  free-stream  Mach  numbers  ranging  from  2.0  to  7.0  in  Ref.  2k. 

The  principal  difference  in  this  model  from  that  of  the  linear  theory  is 

the  elongation  of  Xg,  hence  a  greater  side  force. 

23 

Dugger,  et  al,  at  APL  reexamined  the  same  general  model  shown  in  Fig. 
11  using  Woolard' s  H  results  for  the  forces  due  to  the  shock-flame  system. 
The  model  and  a  typical  pressure  profile  are  shown  in  Fig.  13.  Air  is 
compressed  by  the  bow  shock  and  expanded  by  Prandtl-Meyer  turning  around 
the  knee  before  approaching  the  normal  shock.  The  length  X£  required  for 
the  bounding  streamline  tc  complete  the  turn  is  equal  to  y^  cot  where 
y^  is  the  height  of  the  streamtube  which  will  be  compressed  by  the  normal 
shock  and  captured  by  the  normal  plane  flame  of  height  y^.  If  heat  addition 
is  "supercritical",  the  normal  shock  is  detached  from  the  flame  (Woolard 's 
model),  and  the  flow  will  expand  between  the  shock  and  the  flame,  so  that 
the  flame  height  y^  is  greater  than  y In  the  special  case  of  critical 
(Chapman-Jouguet ,  "lean  limit")  heat  addition,  the  distance  x^  vanishes 
and  y4  =  y2- 
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Fig.  13  SKETCH  OF  WEDGE-TYPE  EXTERNAL  RAMJET  WITH  PLANAR  HEAT  ADDITION 
BEHIND  NORMAL  SHOCK,  AND  TYPICAL  SURFACE  PRESSURE  PROFILE  FOR 
6,  =  62  =  5°,  M0=  5.4,  M2  =  6,  y0  =  Y2  ~  Y4=  '-4- 
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Results  for  flame  height  to  body  thickness  y^/r  and  range  parameter 
9.  are  shewn  in  Fig-  14.  The  triangular  points  are  based  on  the  node!  of 
Fig-  13  using  Woolard’s  approximate  solution  in  the  expansion  region.  The 
solid  triangles  show  that  supercritical  heat  addition  was  required  for 
6,  =  5  >  6°  at  Mrt  =  5  and  for  M-  >  8  with  6.  =  o  =  5°.  These  curves  show 

1  2  —  0  0—  i  2 

that  at  a  flight  Mach  number  of  5,  the  required  flame  height  increases  and 
the  range  capability  of  the  vehicle  decreases  as  wedge  angles  are  increased. 
The  range  falls  drastically  for  wedge  angles  greater  than  6°.  For  engines 
with  6|  =  &2  =  5°,  flace  height  decreases  and  range  increases  as  flight 
Mach  number  is  increased  to  5  or  6,  but  reverse  trends  are  in  evidence  by 
Mach  8- 

The  circular  points  in  this  figure  are  from  the  analysis  of  Fig.  11. 
Both  methods  predict  about  the  same  flame  height  for  zero  net  thrust  but  the 
latter  study  predicts  lower  cruise  range  efficiency  at  higher  Mach  numbers. 
This  is  due  to  the  fact  as  increases,  decreases  and  the  flame  is 
shifted  aft  toward  the  trailing  edge  and  part  of  the  pressure  field  due  to 
combustion  is  lost  in  the  wake-  Likewise,  as  the  wedge  angles  increase  the 
model  becomes  stubbier  resulting  in  the  same  effect-  Clearly  this  would 
be  a  poor  method  of  adding  heat  at  high  Mach  number;  however,  the  problem 
could  be  relieved  simply  by  making  6^  <  5^  within  practical  limits,  thus 
shifting  the  flame  forward.  Also  shown  as  dashed  lines  on  the  figure  are 
results  from  the  constant-pressure  combustion  analysis  from  the  same 
reference,  which  is  discussed  next. 
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2.3  Constant-Pressure  Heat  Addition 

Although  the  planar  flame  models  lead  to  significant  simplification 

in  the  analysis  of  external  burning  and  give  what  appears  to  be  reasonable 

estimates  of  the  integrated  pressure  forces,  they  do  not,  in  general, predict 

the  experimental  pressure  distributions  obtaii  d  in  most  tests.  This  is 

especially  true  for  the  strong  detonation  models  because  the  stationary 

normal  shock  produced  by  heat  addition  alone  has  as  yet  to  be  demonstrated 
* 

experimentally.  Moreover,  as  Mq  increases,  the  losses  across  a  normal 
shock  become  excessive  and  contribute  to  a  decrement  in  performance.  For 
these  reasons,  and  because  considerable  data  had  shown  continuous  heat 
release  at  near  constant  pressure,  the  following  analysis  was  generated. 

Figure  15  shows  the  general  representation  of  the  analytical  treat¬ 
ment  for  continuous  constant-pressure  heat  release.  In  Fig.  15(A)  the  heat 
addition  is  large  enough  to  support  an  oblique  shock  (called  the  flame  shock) 
and  in  Fig.  15(B)  the  heat  addition  is  only  of  sufficient  strength  to  reduce 
the  Prandtl-Meyer  expansion  in  the  flow  around  the  knee.  The  particular  case 

it  26 

The  stabilized  detonation  experiments  conducted  by  Nicolls  used  an 

under-expanded  supersonic  nozzle  with  premixed  fuel  to  produce  the  normal 

,  ,  j  27  28 

shock  and  Gross  and  Rhodes  and  Chriso  generated  a  Mach-reflected  shock 

using  wedges  in  the  side  of  a  two-dimensional  supersonic  tunnel. 
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A)  THEORETICAL  MODEL  FOR  HEAT  ADDITION  LARGE  ENOUGH 
TO  5UPP0RT  AN  OBLIQUE  SHOCK. 


B)  THEORETICAL  MODEL  FOR  SMALL  HEAT  ADDITION 

Fig.  15  ANALYSIS  OF  WEDGE-TYPE  ERJ  WITH  CONSTANT 
PRESSURE  HEAT  ADDITION 
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of  neither  shock  nor  expansion  corresponds  to  a  constant  pressure  field 
on  the  entire  underneath  surface  of  the  model. 


The  analysis  is  based  on  the  following  assumptions: 

(1)  The  air  flowing  around  the  model  outside  the  heated  zone  has 
properties  defined  in  Refs.  29  and  30  ;  the  flow  is  two-dimensional  and 
inviscid.  The  analysis  is  handled  by  stations  in  the  flow  so  that  analysis 
is  actually  one-dimensional. 

(2)  No  heat  or  mass  is  transferred  across  the  interface  between 
the  heated  zone  and  the  adjacent  air  flow. 

(3)  A  given  streamtube  of  air  adjacent  to  the  model  and  within 

the  bow  shock  at  the  model  knee  receives  all  the  fuel.  The  fuel  is  injected 
just  behind  the  flame  shock  in  the  large  heat  addition  case  or  just  behind 
the  rear-running  Mach  line  for  the  Prandtl-Meyer  turn  case  at  condition 


u2’  T2’  p2  (with  Mach  angle  $2)'  *n  eit^er  case,  the  static  pressure  and 
velocity  component  u^  in  the  emergent  fuel  jet  are  matched  to  the  air 
static  pressure  p2  and  air  velocity  u2.  Hence,  per  the  one-dimensional 
treatment  the  following  relations  apply: 


P4  5 

P2 

(30) 

Mass : 

*4" 

w2  +  *f 

(31) 

where 

w2  = 

Wa  =  8P2U2A2/R2T2 

(32) 

Momentum: 

Va 

=  w2u2  +  wfu^  =  (w2  +  wf)  u2 

(33) 

Hence 

u4  = 

U2 

(34) 

and 

*4  = 

gf  2'4^R4^4  ”f  +  8P2U2^2^R2^2 

(35) 

It  is  assumed  that  for  this  analysis  the  maximum  amount  of  air  which  can 


be  used  in  the  combustion  is  represented  by  the  full  amount  captured  by 
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by  the  bow  shock  at  the  model  knee.  This  will  be  referred  to  as  the 
"bow-shock  limit",  for  which 

yl  =  ylmax  =  *1  tan  (91  '  V  =  (Vd)/f  Piul  <36> 

where  d  is  the  width  of  the  two-dimensional  engine,  and  9^  is  the  bow  shock 
angle,  and  x^  =  r/sin  6^. 

(4)  The  side  boundaries  of  the  heated  zone  are  the  model  and  the 

hot-cold  interface.  The  exit  area  (represented  by  y^)  is  in  a  plane 

which  forms  equal  angles  9  with  the  model  surface  and  the  interface.  The 
upstream  boundary  of  the  heated  zone  is  either  the  flame  shock  or  the  last 
Mach  wave  in  the  Prandtl-Meyer  expansion.  The  effective  cross  sectional 
area  normal  to  the  velocity  vector  u?,  at  the  start  of  the  heat  addition  is 
found  for  either  case  by  extrapolating  the  interface  plane  back  to  the  model 
knee  and  is  represented  by  y^*  In  either  case,  constant  pressure  heat 
release  if  '-sumed  to  supply  the  expansion  (by  tempeiature  increase)  to 
bring  the  flow  to  the  area  represented  by  .  The  relationship  between  y^ 
and  y^  is: 

y^  =  y  /con  a  +  2  x^  sin  a  (37) 

(5)  Heat  release  is  completed  at  station  4  with  100%  combustion 
efficiency.  Properties  of  combustion  products  of  kerosene-air  and  hydrogen- 
air  systems  are  taken  from  Refs.  31  and  32  respectively. 

With  the  above  assumptions  (since  velocity  is  constant)  the  energy 
equation  may  be  based  simply  on  static  enthalpies: 

"f2  hf2  '  “4h4  ‘  “2h2  (38) 
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Equations  (32X  (35),  (37)  and  (38)  can  be  combined  to  give  the  fuel  flow 
rate  per  unit  engine  width: 


w^/d  = 


2g  P2U2x2f  h4  sin  a 


T4hff  +  h2T4 


h4T2 


(39) 


^  (cos  a)  ny 

Solution  procedure  is  as  follows.  Conditions  at  station  2  are 

determined  for  various  6^'s  from  two-dimensional  oblique  shock  equations. 

Ther  for  selected  values  of  f,  Eq.  (38)  is  used  to  determine  h.  at  p,  =  p  . 

4  r4  r2 

Other  properties  at  station  4  can  now  be  found  from  combustion  tables.  For 

each  combination  of  f  and  6^  w  is  calculated  from  Eq.  39.  Equation  (36)  is 

then  tested  for  y  <  y  and  limits  on  6  are  established. 

1  ^max  3 

Lift  and  thrust  forces  per  unit  width  are: 

L/d  =  p^  cos  (6^  +  a)  +  p2x2  cos  (§2  -  a)  -  p  c  cos  a 

F/d  =  -  Plxx  sin  (^  +  or)  +  p2x2  sin  (62  -  a)  +  p  c  sin  <y 

For  the  special  case  0=0  and  =  $2>  ^ese  equations  are  simply: 

L/d  = 


h  +  i>2 

' 

2 

"  Po 

F/d  =  (p2  -  pj)  t 


(40) 


//  1  \ 
VHiy 


List  and  thrust  coefficients  are  defined: 


CL  5  Ut<I0  ‘  L/AF  % 

(42) 

Cp  *  F/dTq0  =  F/Ap  ,q 

(43) 

It  is  assumed  that  to  a  first  approximation  the  friction  drag  is  the  same 
with  or  without  burning,  so  that 
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C_  =  (friction  drag  without  combustion) /A,  q~ 

Dj  I  U 

and  the  coefficient  of  net  thrust  is 

CT  =  CF  CDf 

A  cruise  range  parameter  (in  nautical  miles  units)  is  defined: 

2 


ft  h  L  Uq/6076  wf 


1  ■  w 


(44) 


(45) 


(46) 


Results  of  calculations  from  the  above  theory  are  given  in  Fig.  16 
for  kerosene  (C^I^)  fuel,  and  the  effects  of  substituting  hydrogen  fuel 
are  stun-n  in  Fig.  17.  All  curves  presented  are  for  zero  angle-of-attack 
(or  =  0),  because  preliminary  calculations  showed  that  the  optimum  O’  would 
always  be  zero  for  zero  gross  thrust  (C_  =  0)  and  would  always  be  near  zero 
(0°  <  ct  <  3°)  for  small  net  thrusts  (C^,  <  0.1).  Figures  16(A)  and  (B) 
illustrate  effects  of  engine  wedge  angles  at  zero  gross  thrust  and  effects 
of  required  thrust  level  for  a  given  configuration,  respectively,  for 
stoichiometric  combustion  of  kerosene  (ER  =  1).  Figures  16(C)  and  (D)  then 
illustrate  equivalence  ratio  effects,  which  in  turn  determine  the  optimum  6^ 
for  given  thrust  requirement  and  flight  conditions,  as  shown  in  Fig.  16(E), 


A 

The  Breguet  cruise  range  would  be  found  by  multiplying  ft  by  the  natural 
logarithm  of  (gross  weight/burnout  weight).  This  logarithmic  term  is  unity 
when  the  fuel  load  is  63.2%  of  gross  weight,  in  which  case  ft  becomes  the 
range  in  nautical  miles.  The  w^  must  provide  thrust  equal  to  vehicle  drag, 


and  L  must  equal  vehicle  weight.  The 
centrifugal  lift. 


1  -  <»„/»,>* 


factor  corrects  for 
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A)  EFFECT  OF  WEDGE  ANGLES  ON  LIFT  COEFFICIENT  AND  RANGE 
PARAMETER  «  =  0°,  CF  =  0,  q„  =  350  PSF,  KEROSENE,  ER  =  1. 


PSF 
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B)  EFFECT  OF  REQUIRED  THRUST  COEFFICIENT  AND  DYNAMIC 
PERSSURE  ON  RANGE  PARAMETER  M  =  5,  a  =  0°,  S,  =  82  =  5°, 
KEROSENE.  ER  =  I. 


Fis.  16  ANALYSIS  OF  WEDGE-TYPE  ERJ  WITH  CONSTANT  PRESSURE 
HEAT  ADDITION  USING  KEROSENE. 
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C)  yj/yjmox  AS  FUNCTION  OF  ER  AND  AFT  WEDGE  ANGLE  S2 
FOR  ERJ  AT  M0=  9.0  AND  q0=  1000  PSF. 


0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 

EQUIVALENCE  RATIO,  ER 


D>  RANGE  FACTOR  AS  FUNCTION  OF  ER  AND  M0  FOR 
ERJ  WITH  CF  =  0  AND  q0  =  1000  PSF. 


Fig.  16  (Continu«d) 

ANALYSIS  OF  WEDGE-TYPE  ERJ  WITH  CONSTANT  PRESSURE 
HEAT  ADDITION  USING  KEROSENE. 
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(F)  COMPARISON  OF  ERJi  WITH  VARIOUS  THRUST  REQUIREMENTS  TO  CRJ 
VEHICLES  WITH  OVER-ALL  L/D  =  6.  ALL  ENGINES  CRUISE  WITH  KERO 
SENE  FUEL  AT  OPTIMUM  ER. 


Fig.  16  (Concluded) 

ANALYSIS  OF  WEDGE-TYPE  ERJ  WITH  CONSTANT  PRESSURE 
HEAT  ADDITION  USING  KEROSENE. 
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and  hence  the  optimized  or  on-design  ft  vs  Mr  characteristics  as  shown  in 
Fig.  16  (F).  Notable  conclusions  and  explanatory  r.-.-aarks  from  these  figures 
are: 

(1)  From  Fig.  16(A)  it  is  seen  that  cruise  range  parameter  for  a 
frictionless  cruise  vehicle  (C  r-  0)  decreases  slightly  as  aft  wedge  angle 

is  increased  from  5°  to  30°  for  a  given  fore  wedge  angle  6^. 

(2)  For  a  given  6^,  the  cruise  range  parameter  increases  as  6^  is 
decreased,  because  the  bow  shock  compression  is  decreased,  the  wave  drag 
is  correspondingly  reduced,  a  small  air  mass  flow  is  ’captured"  by  the 
bow  shock,  and  the  net  result  is  that  less  fuel  is  required  to  produce 

C„  =  0  by  stoichiometric  combustion.  The  total  lift  force  L  decreases,  but 
the  requited  fuel  rate  decreases  faster,  so  that  L/w^  and  ft  increase. 

(3)  Lift  coefficient  referred  to  frontal  area  per  Eq.  (42)  decreases 

as  increased,  because  the  pressure  due  to  burning  is  felt  on  a  smaller 

planform  area  and  the  total  lift  force  for  an  engine  of  fixed  frontal  area  is 

reduced.  (Lift  coefficient  referred  to  planform  area  would  be  the  same  for 

all  6  when  C„  =  0  and  a  =  0°  since  in  this  case  the  lift  is  proportional 
4  F 

to  planform  area). 

(4)  It  is  concluded  from  the  above  and  from  other  similar  calcula¬ 
tions  that  maximum  ranges  will  be  obtained  with  slender  configurations  having 
the  smallest  practical  6^  which  probably  will  be  near  5°.  Remaining  figures 
are  for  6^  =  5°  and  >5°,  depending  upon  required  thrust  level. 

(5)  The  effect  of  required  thrust  level  is  illustrated  by  the  plot 

of  ft  versus  C  in  part  B  for  6,  =  6  =  5°.  It  is  seen  that  ft  decreases 

F  1  4 
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rapidly  as  required  thrust  is  increased,  even  to  tht  amount  required  to 

•k 

overcome  friction  drag.  This  is  probably  the  most  significant  general 

result  of  the  theory,  because  it  suggests  that  ERJ*k  probably  will  not 

be  attractive  for  high  thrust  or  accelerating  missions.  It  is  discussed 

further  in  connection  with  Figs.  16(E),  16(F)  and  17. 

(6)  Figure  16(B)  also  shows  that  the  effect  of  flight  dynamic 

pressure  q^  is  small  per  this  analysis.  (As  discussed  elsewhere,  experimental 
factors,  such  as  boundary  layer  effects  and  non-  ideal  fuel  distribution,  may 

lead  to  much  stronger  real  effects  of  q^  or  flight  altitude). 

(7)  Figure  16(D)  shows  that  for  a  given  configuration  at  given  flight 

conditions,  ft  increases  as  ER  is  decreased.  However,  the  extent  to  which  ER 

can  be  decreased  is  limited  in  the  present  analysis  by  the  criterion  y^  < 

yl  max  per  E°*'  ^)*  Figure  16(C)  shows  how  these  "bow  shock  limits"  are 

established  for  various  thrust  levels  and  at  Mq  =  9.  It  is  seen  that 

ER  decreases  as  y,/y.  increases;  the  optimum  or  bow-shock-limited  value, 

which  represents  the  minimum  ER  and  hence  the  maximum  ft  by  these  rules,  is 

that  corresponding  to  y^/y^  ^  ■=  1.0.  Typical  values  of  this  optimum  ER 

for  various  conditions  are  given  in  Table  II. 

*  The  friction  drag  coefficients  calculated  for  part  (D)  and  subsequent 
parts  are  based  on  a  50-ft  chord;  for  similar  geometries,  longer  engines  would 
have  smaller  's  and  vice  versa. 
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TABLE  II 

Bow-Shock-Limited  Values  of  Kerosene  Equivalence  Ratio**  for  Various  Flight 
Mach  Numbers,  Thrust  Levels,  and  S^s  for  ERj's  of  50-ft  Chord  at  =  1000  psf. 
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As  one  would  expect,  Table  II  shows  that  must  be  increased  as  required 

thrust  is  increased.  For  given  thrust  level  and  engine  length,  the  minimum 

ER  decreases  as  6„  is  increased,  because  frontal  area  and  y,  increase  as 

the  maximum  thickness  point  is  shifted  toward  the  rear  of  a  fixed-length 

engine.  The  latter  effect  is  more  pronounced  at  positive  net  thrust. 

(8)  Figure  16(E)  shows  that  an  optimum  6^  (with  its  corresponding 

minimum  ER)  will  exist  for  any  thrust  level  at  given  flight  conditions. 

For  C„  =  0,  6„  opt  ~  6  .  For  higher  thrusts,  6„  _  >  6,.  Such  optimum 

F  2  r  —  1  2  oot  1  r 


ER's  less  than  0.1  or  greater  than  1.0  are  considered  unrealistic.  In 
fact,  ER’s  <  0.3  may  not  be  realistic  for  kerosene,  but  there  might  be  some 
advantage  in  going  slightly  rich,  say  to  ER  ~  1.2. 
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o^'s  are  used  to  plot  maximum  ft  versus  Mg  for  various  thrust  levels  in 
Fig.  16(F).  The  drastic  deterioration  in  ft  as  required  thrust  is  increased 
is  clearly  seen  in  these  two  figures.  For  example,  from  Fig.  16(F)  it  is 
seen  that: 

(a)  ft  is  8300  n.m.  at  Mg  =  15  (and  still  rising)  if  Cp  =  0, 
in  which  case  6^  =  6^  =  5°  is  optimum 

(L)  ft  is  reduced  to  520u  n.m.  at  Mn  =  15  for  =  0  for  50-ft 

engines  with  5.  or  to  4500  n.m.  if  5  is  held  at  5°,  and 
°  2  opt  2 

(c)  if  required  is  raised  to  0.1  to  accommodate  external 
drag  of  other  vehicle  components  or  maneuver  "noise",  ft  decreases  continuously 
from  2000  n.m.  at  Mg  =  5  to  900  n.m.  at  Mg  =  12. 

(9)  ERJ's  should  be  superior  to  conventional  ramjets  with  subsonic 
internal  combustion  (CSJ's)  in  „ruise  applications  at  flight  Mach  numbers 
above  9  or  so,  depending  on  the  assumptions  made  for  the  two  engines.  Two 
CRJ  curves  are  given  in  Fig.  16(F).  Both  are  based  on  engines  with  inlet, 
combustion,  and  nozzle  efficiencies  of  0.92,  0.95  and  0.96,  respectively, 
and  on  vehicles  with  over-all  L/D's  of  6.  The  more  optimistic  CRJ  curve  is 
based  on  chemical  equilibrium  in  the  exhaust  flow,  whereas  the  poorer  curve 
is  based  on  estimated  "actual"  exhaust  flow  properties.  The  latter  flow 
properties  are  intermediate  between  equilibrium  and  frozen  exhaust  flow 
properties  and  are  estimated  as  in  Ref.  33,  The  most  reasonable  comparison 
of  ERJ  and  CRJ  from  this  figure  is  probably  that  between  the  ERJ  curve  for 
C.J,  =  0  and  6^  t  and  the  "actual"  flow  CRJ  curve.  On  this  basis,  the  ERJ 
is  superior  for  Mg  >~8.5.  Longer  ERJ's  would  have  lower  friction  drag 
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coefficients  and  slightly  improved  performance  Or  vice  versa. 

It  should  be  pointed  out  that  the  ERJ  curves  are  (as  done  here) 
properly  based  on  equilibrium  flow  properties,  because  it  has  been  assumed 
that  heat  is  added  continuously  to  the  trailing  edge  of  the  model  and  no 
credit  has  been  taken  for  subsequent  expansion.  (A  similar  assumption  for 
the  CRJ  would  permit  only  a  convergent  nozzle,  and  performance  would  be 
much  poorer  than  the  "actual"  flow  curve  presented  for  the  CEJ.)  Further¬ 
more,  the  analytical  model  of  Fig.  15(A),  is  based  on  the  assumption  that 
heat  is  added  all  the  way  across  a  one-dimensional  flow  to  the  exit  plane 
A^.  Actually,  the  heat  added  downstream  of  the  last  forward-running  Mach 

•k 

wave  striking  the  trailing  edge  cannot  affect  the  pressure  on  the  model 
surface,  hence  it  is  wasted.  A  controlled  fuel  distribution  and  heat  re¬ 
lease  program  might  therefore  decrease  the  specific  fuel  consumption  and 
significantly  improve  the  relative  ERJ  performance  at  all  thrust  levels 
illustrated. 

(10)  Both  CRJ's  and  ERJ's  have  much  greater  cruise  ranges  when 
hydrogen  fuel  is  used  instead  ;f  kerosene,  as  shown  in  Fig.  17.  For  a 
given  fuel  mass  fraction  .-.nd  cruise  Mach  number,  the  range  is  increased 
by  approximately  a  fuctor  of  2.8.  Thus,  an  ERJ  with  an  initial  hydrogen 
load  equal  to  63 %  of  its  gross  weight  could  reach  10,500  n.rn.  at  Mach  12 
on  a  zero  net  thrust  cruise. 


*  This  refers  to  the  left-running  characteristic  of  the  burned  gas  flow  as  it 
passes  the  trailing  edge  of  the  model.  The  pressure  field  in  the  supersonic 
flew  downstream  of  this  Mach  wave  cannot  be  transmitted  upstream. 
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IWc-  additional  Refinements  in  the  constant  pressure  heat  addition 
model  were  made  by  Billig  in  Ref.  34.  These  changes  are  the  inclusion  of 
a  separated  z on?  rn  the  region  of  the  fuel  ports  and  an  approximation  for 
the  pressure  decay  downstream  of  the  heat  release  zone  as  shown  in  Fig.  13, 
and  they  are  more  apparent  when  the  analysts  is  applied  to  a  case  of  heat 
addition  adjacent  to  a  flat  plate  but  are  applicable  to  the  wedge  model. 

Fuel  is  injected  near  the  leading  edge,  causing  a  separated-flow 

region  and  a  corresponding  oblique  injection  shock.  The  fuel  penetrates 

to  a  height  and  combustion  begins  in  the  plane  labeled  and  is  completed 

at  plane  Yy  For  additional  simplicity  the  heat  release  is  assumed  to 

proceed  at  a  rate  to  maintain  the  pressure  level  attained  by  the  injection 

separation.  Thus,  the  injection  shock  is  represented  by  a  straight  line 

until  point  C  is  reached,  where  expansion  commences.  At  point  D  in  the  Y^ 

plane,  expansion  begins  in  the  combustion  gas,  as  represented  by  the  "leading 

Mach  lines",  DE.  Pressure  is  therefore  assumed  to  be  constant  throughout 

the  region  ABDE,  so  that  the  surface  pressure  rise  over  the  length  from 

A  to  E  is  (p£  -  Pq).  To  avoid  the  complex  characteristic  solution  of  the 

flow  field  of  the  combustion  gases,  the  surface  pressure  in  the  expansion 

region  (EF)  is  assumed  to  be  the  arithmetic  mean  of  combustion  and  free 

stream  pressures,  i.e.,  p__  =  (p„  -  p_)/2.  Equilibrium  combustion  of  tri- 

br  2  U 

ethyl  aluminum  in  air  with  1007»  combustion  efficiency  using  the  procedure 

35 

developed  by  Browne  and  Williams  was  carried  out  on  a  high-speed  computer. 
No  heat  or  mass  transfers  across  the  hot -cold  interface. 

The  normal  force  specific  impulse  is  defined  as  the  total  force  (per 
unit  plate  width)  resulting  from  pressure  above  ambient  divided  by  the  fuel 
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flow  rate  (per  unit  plate  width): 

Fs  (p2  -  PQ)  X4  +  (X5  -  X,)  (p2  -  pq)/2  (p2  -  p0)  _  X4  +  X5 

l™  =  *~£  =  gf  5R  PqVj  =gf  ER  P0u0  2  Y1 

The  norma1  force  coefficient  is  the  total  side  force  divided  by 
the  product  of  the  area  subjected  to  pressure  above  ambient  (this  area 
per  unit  plate  width  is  X,.)  and  free  stream  dynamic  pressure: 


(47) 


C  =JL 

N  V 


(p2  -  P(.)  r (x4/x5)  +  i] 
2  % 


(48) 


Table  III  summarizes  the  computer  results  for  cases  having  a  range 
or  free  stream  Mach  numbers  (Mq),  altitudes  (2),  equivalence  ratios  (ER) , 
and  injection  shock  strengths  (0^,  as  determined  by  6^,  the  air  deflection 
angle  due  to  the  separation  zone).  Cases  labeled  with  asterisks  have  injec¬ 
tion  shock  strengths  and  pressure  ratios  which  correspond  to  separation 

36 

of  a  turbulent  boundary  layer  as  predicted  by  Mager.  Experiments  with 

37,38 

non-reacti.ve  injection  have  resulted  in  reasonable  verification  of 

the  average  pressure  as  estimated  by  Mager.  In  the  subsequent  discussion  of 

experiments  of  triethyl  aluminum  injection  onto  a  flat  place  at  Mach  5 

from  Ref.  34  f  tha  injection  shock  was  somewhat  weaker  and  separation  was 

apparent.  Changes  in  the  geometry  of  the  fuel  injector  may  permit  control 

of  6^,  so  that  all  of  the  cases  listed  can  be  of  interest.  The  cases 

studied  were  limited  co  those  in  which  the  Mach  number  at  the  end  of 

combustion  was  supersonic  (M^  >  1);  however,  subsonic  solutions  are 

possible  and  may  in  fact  be  desirable  at  low  flight  Mach  numbers.  (The 

I..,  and  C„  columns  in  Table  III  will  be  discussed  later.) 
fN  N 
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In  Fig.  19  the  I _ and  C„  dependence  on  M-  is  shown  for  ER's  of 

Zc1  h  0 

1.0  and  0.25  a£  35.0*50  ft  altitude.  For  lean  mixtures,  I is  about  20CC 

fS 

seconds  higher.  For  deflections  equal  to  the  defined  turbulent  boundary 

layer  separation,  6^  -  6^  sep»  c^e  iEpulse  is  essentially  invariant  with 

M_.  Even  for  the  cases  of  constant  6,  =  5°,  variation  of  I_,  with  is 
v  1  za  v 

small,  a  conotonic  decrease  with  increasing  Mq-  The  norcal  force  coefficient, 

C„,  is  mainly  dependent  on  deflection  angle,  increasing  only  slightly  with 

lower  HR.  As  for  cost  aerodynamic  surfaces,  larger  deflection  results 

in  a  greater  force  coefficient. 

Figures  20(A)  and  (B)  show  the  ER  and  deflection  effects  on  I 

more  clearly  for  cases  in  which  all  other  parameters  are  held  constant. 

Figure  20(C)  shows  the  altitude  effect  at  =  5,  6,  =  6,  .  At  constant 

0  1  1  sep 

Mq  (but  varying  u^) ,  temperature  rather  than  pressure  has  the  dominant 
effect  on  1^,  as  shown  by  the  rapid  rise  in  1^  from  sea  level  to  the 
tropopause  (Z  -  36K  ft)  and  correspondingly  slight  decline  in  the  constant- 
temperature,  decreasing-pressure  region.  Of  course,  if  a  comparison  were 
made  on  a  constant  velocity  basis  (say,  =  5000  ft/sec),  the  effect  of 
altitude  would  be  smaller  (but  deflection  angles  would  have  to  be  adjusted 
for  Mq  variation).  ICAO  Standard  Day  air  properties  were  used  throughout. 

The  results  of  the  theoretical  study  are  dimensionless  with  regard 
to  the  geometric  variables.  For  flight  models  certain  constraints  will 
bound  the  range  of  practical  geometries.  Two  principal  limitations  will 
be  a  minimum  practical  combustor  length  for  efficient  combustion  and  a 
maximum  height  of  fuel  jet  penetration.  The  latter  effect  is  more  pronounced 
for  low  ER  ratios  since  at  the  same  combustion  pressure  low  ER  requires 
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Fig.  19  NORMAL  FORCE  COEFFICIENTS  AND  SPECIFIC  IMPULSES  FOR 
A  FLAT  PLATE  WITH  CONSTANT  PRESSURE  COMBUSTION  OF 
TEA  FUEL 
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Fig.  20  NORMAL  FORCE  SPECIFIC  IMPULSE  AS  FUNCTION 
OF  EQUIVALENCE  RATIO,  DEFLECTION  ANGLE 
AND  ALTITUDE  FOR  FLAT  PLATE  AT  MACH  5 
WITH  CONSTANT  PRESSURE  COMBUSTION  OF  TEA 
FUEL 
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relatively  greater  jet  penetration  to  obtain  the  same  length  of  the 
constant  pressure  field  A-E.  This  effect  is  dramatically  shown  in 
Table  III.  The  last  two  columns  in  Table  III,  I£„  and  C„,  are  the 
impulse  and  force  coefficients  which  would  result  if  the  fuel  penetration 
depths  were  limited  to  the  theoretical  value  for  stoichiometric  fuel 
injection,  i.e.,  Y.  =  Y,  »  _  , .  In  effect,  this  constraint  forces 

the  expansion  to  be  completed  before  the  end  of  the  plate  is  reached 
when  ER  <  1.  The  flow  model  would  then  include  the  distance  X£  (Fig. 18) 

D 

with  zero  pressure  coefficient  (C„  =  0),  hence  the  over-all  C„  would 

>  6  N 
decline,  and  the  1^  **  1^  produced  on  the  constrained  length  X,.  would 

decrease.  The  net  result  reverses  the  trend  of  better  performance  for 
lower  ER  and  instead  shows  the  desirability  of  rich  operation  if  pene¬ 
tration  rather  than  length  is  controlling. 

The  above  theoretical  treatments  can  be  applied  directly  to 
planar  surfaces  at  angles-of-attadc  (i.e.,  other  than  zero  incidence 
to  the  free  stream).  It  can  also  be  readily  modified  to  handle  single¬ 
port  rather  than  multiport  or  iine-source  injection  by  using  half-conical 
flow  surfaces  and  deflection  surfaces  rather  than  two-dimensional  ones. 

39 

In  closing  this  summary  of  theoretical  studies,  the  work  or  Wald 

40-43 

of  United  Aircraft  and  Yen  and  McCloy  of  General  Dynamics  should  be 

mentioned.  In  particular,  fundamental  theoretical  studies  by  McCloy 

(Department  of  Aeronautical  and  Astronautical  Engineering,  University  of 

Illinois  and  consultant  to  APL  in  1960-61)  led  to  the  flow  model  concepts 

used  in  the  preceding  constant  pressure  analysis.  Noteworthy  too,  are  the 
44 

studies  by  Marino  who  apparently  adapted  McCloy' s  model  to  the  flat 
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plata  case  and  studied  both  the  simple  two-dimensional  (wedge  heat  addition) 
and  a  conical  heat  addition  (i.e.,  constant  pressure  combustion  from  a  point 
source)  and  concludes  that  the  former  is  more  efficient.  He  presents  per¬ 
formance  charts  for  =  2  to  20  and  came  to  similar  conclusions  to  those 
showed  in  Ref.  36. 

3.  LITERATURE  REVIEW  OF  EXPERIMENTAL 
STUDIES  OF  EXTERNAL  BURKING 

The  first  experiments  related  to  external  combustion  were  concerned 

with  the  reduction  of  base  drag  of  projectiles  by  combustion  in  the  wake. 

45 

Baker,  et  al  injected  hydrogen  into  the  base  of  a  2-1/4  inch  diameter 

cone-cylinder  placed  in  a  Mach  1.6  free  jet  and  obtained  base  drag  reduction 

of  60-75%.  Subsequent  additional  analysis  of  the  data  suggests  that  the 

specific  impulses  based  on  drag  reduction  were  1200  to  4000  secs  and  the 

combustion  efficiencies  were  80-90%.  Similar  tests  were  made  by  Scanland 

4-6 

and  Hebrank  who  burned  a  solid  pyrotechnic  composition  in  the  base  of 
40  mm  projectiles.  The  base  drag  reduction  at  Mq  =  1.85  was  657«  which 

corresponded  to  a  total  drag  reduction  of  19%.  In  both  of  these  tests,  the 

burning  was  external  but  was  probably  confined  to  the  subsonic  portions  of 
the  wake. 

Davis,  at  Experiment  Incorporated,  made  a  series  of  tests  with  hydrogen 
injection  from  a  flat  plate  at  =  1.7.  Typical  pressure  rises  due  to  com¬ 
bustion  are  shown  in  Fig.  21  for  hydrogen  injection  from  the  base  of  a  rear¬ 
ward  facing  step.  In  nearly  all  tests  not  having  a  high  drag  f lameholder 

immersed  in  the  external  air  system  the  addition  of  oxygen  was  necessary  to 
establish  burning.  Even  with  oxygen  addition  it  is  apparent  (e.g.,  burning 
in  Run  a  but  not  in  Run  b) 
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Fig.  21  EXTERNAL  COMBUSTION  MODEL  TESTED  BY  DAVIS  ET  AL.  WITH  TYPICAL 
STATIC  PRESSURE  PROFILES 
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that  the  combustion  is  near  a  stability  limit.  The  difficulties  encountered 

in  getting  hydrogen  to  burn  in  these  tests  at  1  atm  and  about  5rO°R  static 

are  not  at  all  surprising.  In  fact,  on  the  basis  of  recent  ignition  delay 
47 

experiments,  it  is  doubtful  that  ignition  of  hydrogen  could  have  occurred 
unless  a  significant  separated  zone  with  recirculation  existed  downstream 
of  the  step.  A  large  number  of  injector  and  flame-holder  configurations 
were  tried,  some  with  supplemental  oxygen  and  most  with  spark  ignition  source. 
Although  most  of  the  data  are  difficult  to  interpret  due  to  the  presence  of 
flame-holder  shock  waves,  a  general  conclusion  can  be  drawn  that  when  com¬ 
bustion  occurred  the  "volume  source"  effect  essentially  compensated  for  the 
expansion  due  to  the  step,  and  thus  represented  a  net  side  force  increase 
due  to  combustion. 

On  the  basis  of  some  of  the  aforementioned  theoretical  studies, 

48-53 

Dorsch,  et  al  reasoned  that  to  obtain  the  desired  performance  from 

external  burning  it  would  be  necessary  to  stabilize  combustion  in  the  super¬ 
sonic  portion  of  the  flow  field  without  the  attendant  drag  losses  of  flame- 
holders,  as  were  necessary  in  many  of  Smith  and  Davis'  tests.  However, 
since  the  conditions  for  combustion  were  unfavorable;  i.e.,  near  free-stream 
ambient  pressure  and  temperature  for  the  flight  condition  with  residence 
times  of  less  than  100  ps,they  felt  that  conventional  fuels  were  not  practical. 
For  this  reason,  they  used  very  reactive  fuels,  principally  aluminum  borohydride, 
in  their  tests  and  it  is  notable  that  the  only  successful  external  burning 
experiments  made  since  that  time  have  been  with  reactive  fuels.  In  Ref.  48} 
aluminum  borohydride  was  injected  from  the  walls  of  a  3.84-inch  by  10  inch 
supersonic  wind  tunnel  into  a  Mach  1.5  to  4.0  airstream.  Since  the  maximum 

tunnel  plenum  conditions  were  limited  to  540°  to  570°R  and  1-1.5  atm,  the 
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static  temperature  and  pressure  in  the  test  section  decreased  with  increas¬ 
ing  Mach  number,  dropping  from  386°R  and  0.4  atm  at  Mach  1.5  to  133°R  and 
0.01  atm  at  Mach  4.0.  In  spite  of  the  high  altitude  pressure  simulation 
(7lKft  @  Mq  =  3)  and  unrealistically  low  temperatures,  (~200°R  low  @  Mq  =  3) 
steady  combustion  was  demonstrated  at  Mq  =  1.5,  2  and  3.  Spontaneous  ignition 
occurred  in  the  Mach  1.5  and  2.0  tests,  but  at  Mach  3.0,  a  downstream  spark 
source  was  needed  which,  after  ignition  and  flashback  to  the  injector,  could 
be  terminated.  In  Ref .49  ,  wall  static  pressures  were  measured  in  the  same 
facility  at  Mq  -  2  and  3.  Pressure  rises  of  from  20  to  40  percent  above  the 
no  injection  values  were  obtained  which  had  an  increasing  trend  with  increas¬ 
ing  equivalence  ratio  and  slightly  higher  ratios  at  Mq  =  3  than  at  MQ  =  2. 

These  results  prompted  an  investigation  of  a  flat-plate  model  in  a  larger 

50 

Mach  2.46,  1  ft  by  1  ft  tunnel.  In  this  facility,  the  nominal  static  pres¬ 
sure  was  0.1  atm  and  the  static  temperature  was  250°R  in  the  absence  of 
combustion.  A  schematic  illustration  of  a  typical  flame  shape,  shock  wave 
pattern  and  chordwise  pressure  distribution  are  shown  in  Fig.  22(A).  The 
model  consisted  of  a  short  basic  flat  plate  with  a  12-in  span  and  13-in 
chord  plus  an  extension  plate  with  an  additional  12-in  of  chord.  As  sus¬ 
pected  by  the  authors  in  their  analysis  of  the  data  and  later  verified  in 
a  larger  (10-ft  x  10-ft)  tunnel  test  of  a  similar  model, ^  the  secondary 
pressure  rise  downstream  was  due  to  tunnel  effects.  A  comparison  of  static 
pressure  distributions  from  the  two  facilities  is  shown  in  Fig.  22(B).  Com¬ 
puted  Mach  numbers  based  on  static  and  pitot  pressures  taken  within  the 
luminous  flame  zone  showed  that  the  Mach  number  was  subsonic  through  most 
of  the  flame  zone  and  became  sonic  and  low  supersoni  as  the  hot-cold  inter- 
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Fig.  22  RESULTS  OF  FLAT-PLATE  TESTS  IN  SMALL-  AND  LARGE-SCALE 
TUNNELS  EIY  NASA.  (Ref*.  50  and  51) 
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face  was  approached.  An  analysis  of  this  test  will  be  discussed  later. 

A  body  of  revolution  (a  10.5-inch  parabolic  forebody  followed  by 

a  10.5-inch  cylindrical  afterbody  with  a  1.75-inch-diameter)  was  also 

tested  in  the  1-ft  by  1-ft  Mach  2.46  tunnel.  Fuel  was  injected  through  four 

circumferential  and  equally-spaced,  0.015-inch-diameter  fuel  orifices  which 

were  located  in  the  plane  of  the  forebody-afterbody  junction.  As  was  the 

48-53 

case  in  all  of  the  tests  in  this  series  no  attempts  were  made  to  vary 

the  fuel-flow  rate  or  to  control  or  measure  it  to  a  high  degree  of  accuracy. 

Spark  ignition  was  required.  Combustion  produced  a  near  constant  pressure 

field  on  the  afterbody  with  a  rise  of  about  1  psia  or  about  707„  above  the 

non-burning  condition.  A  corresponding  rise  in  base  pressure  from  0.68 

psia  to  2.5  psia  also  occurred  but  quantitatively  undetermined  wind  tunnel 

effects  may  have  influenced  this  result. 

The  remaining  tests  in  this  series  were  made  with  a  two-dimensional, 

6%  thick,  blunt-base  airfoil,  at  2°  angle-of-attack  in  the  1-ft  by  1-ft 
53 

tunnel.  The  model  had  a  13-inch  chord  and  was  tested  at  Mq  -  2.47  and 
2.96.  The  location  of  th°  heat  addition  region  below  the  wing  was  arbitrarily 
chosen  to  provide  primarily  a  lift  increase  rather  than  a  drag  reduction,  in 
order  to  facilitate  comparison  with  the  theoretical  calculations  of  Ref.  11. 
Figure  23^1)  is  a  combination  open-shutter  and  schlieren  photograph  of  the 
flame  and  associated  shock  wave  system  for  a  test  at  Mq  =  2.47.  In  Fig.  23(B) 
the  typical  pressure  distributions  with  and  without  combustion  for  both  sur¬ 
faces  at  each  Mach  number  tested  are  shown.  Integration  of  the  pressure-area 
distribution  produced  the  forces  and  coefficients  shown  in  Table  IV. 
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A)  COMBINATION  OPEN  SHUTTER  AND  SCHLIEREN  PHOTOGRAPH 
OF  COMBUSTION  AT  M0  =  2.47 


DISTANCE  FROM  LEADING  EDGE,  RATIO  TO  CHORD  LENGTH 
B)  SURFACE  STATIC  PRESSURE  DISTRIBUTIONS 

Fig.  23  RESULTS  FROM  EXTERNAL  BURNING  TESTS  OF  AIRFOIL 
AT  M0  =  2.47  AND  2.96  BY  NASA.  (Ref.  52) 
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TABLE  IV 

o 

Summary  of  Forces  on  Airfoil  at  or  =  2 


Mach 

Condition 

Lift 

c. 

Drag  Force,  lb/ft  span 

L/D 

Force, 

lb/ft 

span 

L 

Upper 

Lower 

Base 

Friction 

Total 

2.47 

Non-burning 

60.7 

0.06 

Burning 

117.7 

0.12 

2.96 

Non-burning 

50.3 

0.05  -1.4 

11.5 

-1.4 

3.2 

11.9 

4.2 

Burning 

113.0 

0.12  -1.4 

15.7 

-2.8 

3.2 

14.7 

7.7 

The  friction  drag  was  based  on  a  skin-friction  coefficient  (on  planform  area)  of 
0.0035  obtained  from  unpublished  wind  tunnel  tests.  The  lift  force  with  combus¬ 
tion  about  doubled  at  both  Mach  numbers.  At  Mq  =  2.96,  which  was  the  only  Mach 
number  that  drag  data  were  analyzed,  the  drag  also  increased  with  combustion. 

This  drag  increase  results  from  an  increased  positive  pressure  coefficient  of 
the  forward  facing  portion  of  the  lower  surface  which  more  than  compensates  for 
the  increased  pressure  on  the  rearward  facing  aft  surface  and  on  the  base. 

With  the  lift  coefficient  defined,  then  the  lift  parameter 

2  A 

C.  (M  -  1)K  is  known,  and  it  is  only  necessary  to  define  an  effective  heater 
0 

deflection  angle  6  ,  and  the  fraction  of  chord,  f,  below  which  the  heat  is  added, 

to  compare  the  experimental  results  with  the  theoretical  analysis  of  section  2. 

Within  the  accuracy  of  the  measurements,  6  was  found  to  be  about  10°  at  both 

n 

Mach  numbers  and  values  of  f  =  0.6  at  Mq  =  2.47  and  f  =  0.8  at  2.96  were 
estimated.  Using  these  values  to  obtain  the 'effective  lift  flow  deflection 
parameter,  a  +  ^  f  6  ,  the  experimentally  determined  lift  parameter  is  compared 
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in  Table  V  with  Che  linear  theory  values  and  the  graphically  determined  values 
for  the  5°  biconvex  airfoil  studied  by  Pinkel,  et  al**  (interpolated  from  points 
shown  in  Fig-  5). 

TABLE  V 

Comparison  of  Measured  and  Theoretical  Lift  Parameters  for  Airfoil 


C'  t-  \  f  6 


Measured  Linear  Theory'  Graphical 


2.47 

5° 

0.28 

0.35 

0.28 

2.96 

6° 

0.34 

0.41 

0.30 

The  significance  of  the  rather  close  correlation  of  theory  with  experiment  is 

that  the  selected  value  of  5  is  indeed  reasonable,  and,  therefore,  a  reasonable 

H 

evaluation  of  the  efficiency  based  on  the  theoretical  methods  for  the  heat  addi¬ 
tion  process  can  be  made.  Equation  (9)  can  be  used  to  find  the  effective  q^ 
which  can  be  ratioed  to  th^  maximum  (complete  combustion)  q,^  available  in  the 
fuel.  In  these  tests,  the  average  wf  was  0.015  ib/sec  which  for  a  heating  value 
of  24,800  Btu/lb  for  aluminum  borohydride  gives  the  efficiencies  listed  in 
Table  VI. 


TABLE  VI 

Heating  Rates  and  Efficiency  for  Airfoil  at  a  =  2° 


Mo 

Ib/sec  q 

2 

.  Btu/sec  ft 

A 

(Actual) 

9t 

2 

Btu/sec  ft 

(Available) 

Efficiency  =  9A/qT 

2.47 

0.015 

311 

575 

0.54 

2.96 

0.015 

226 

431 

0.52 

Thus 

,  about  half 

of  the  effective 

heal 

release  was 

obtained  in  these  tests. 
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Even  though  it  is  evident  in  all  of  these  tests  that  the  heat  addi¬ 
tion  occurs  in  a  zone  which  is  spread  out  in  the  streamwise  direction  and 
therefore  cannot  be  considered  as  a  planar  flame,  it  is  of  interest  to 

examine  the  results  using  a  form  of  the  planar  heat  addition  models  to  see 

54 

whether  any  correlations  can  be  obtained.  Woolard  studied  the  flat  plate 
results  from  Refs.  50  and  51,  using  the  pressure  profiles  from  the  10-ft 
by  10-ft  tests  (Fig.  22(B))  and  the  detailed  in-stream  measurements  from  the 
1-ft  by  1-ft  tests.  Instead  of  the  simpler  single  heated  streamtube  concept 
of  Ref.  18,  he  postulated  the  two-layer  model  shown  in  Fig.  24.  In  the  inner 
layer,  the  Mach  number  just  behind  the  oblique  planar  heat  addition,  M^  is  sub¬ 
sonic.  The  flows  in  regions  3  and  3  are  assumed  to  be  uniform,  parallel  and  nave 
equal  static  pressure.  With  no  further  constraints,  there  are  enough  free 
parameters  to  permit  a  favorable  preferential  matching  of  theoretical  and 
experimental  values.  Those  properties  for  which  an  attempt  was  made  to 
obtain  an  approximate  matching  were  the  peak  pressure  rise,  the  downstream 
pressure  decay,  the  over-all  flame  shape,  and  the  total-temperature  increase 
in  the  inner  layer.  These  results  are  summarized  in  Fig.  25.  The  experi¬ 
mental  flame  shape  was  specified  as  the  luminous-non- luminous  boundary  as 
determined  from  the  sketch  shown  in  Fig.  22(A)  which  apparently  is  an 
unrealistic  description  of  the  method  of  heat  addition.  In  Fig.  25(A)  the 
correlation  between  the  experimental  and  theoretical  pressure  distributions 
and  "planar  flame  shapes"  are  close,  but  the  calculated  temperature  and 
Mach  number  distributions  at  Station  I  of  Fig.  22(A)  shown  in  Fig,  25  (B) 
and  (C)  are  poor,  although  better  than  a  linear  theory  approximation. 
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DISTANCE  PROM  LEADING, EDGE  (inchot) 

A)  PRESSURE  DISTRIBUTIONS  AND  FLAME  SHAPES 


AT,/t,q 


B)  TOTAL  TEMPERATURE  DISTRIBUTIONS  AT  STATION  I 


C)  MACH  NUMBER  DISTRIBUTIONS  AT  STA1  ION  I 


Fig.  25  COMPARISON  OF  FLAT-PLATE  EXPERIMENTAL  RESULTS  WITH 
TWO-LAYER  ANALYSIS  OF  WOOLARD 
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Testing  of  external  burning  ramjets  at  the  Applied  Physics  Laboratory 
v?as  initiated  in  1958.  The  first  phase  of  the  program  consisted  of  small- 
scale  wedge  model  tests  in  a  6-in.  by  7-in.  high-temperature,  Mach  5.0 
propulsion  tunnel  with  100K  ft  altitude  simulation.  The  principal  objective 
in  the  first  phase  of  the  program  was  to  demonstrate  useful  thrust  on  a 
configuration  similar  to  the  type  shown  in  Fig.  2(C).  In  the  second  phase, 
a  scaled  version  of  one  of  the  smaller  models  was  tested  in  a  32-in. -diameter 
free  jet  at  the  Ordnance  Aerophysics  Laboratory  (OAL)  at  Mach  5.0  with  some¬ 
what  less  than  simulated  static  temperature  (260°R  vs  390°R)  and  with  pres¬ 
sure  altitude  simulation  from  82K  ft  to  105K  ft.  In  the  last  phase  of 
the  program,  extensive  in-stream  measurements  were  made  on  a  flat  plate  model 

tested  in  a  10-in. -diameter  Mach  5.0  tunnel  with  true  simulation  of  flight 

34 

at  66K  ft.  Aluminum  alkyl  fuels,  which  reacted  spontaneously  without 

supplemental  ignition  sources,  were  used  in  the  tests.  Fuel  handling,  metering 

34 

and  heating  systems  were  developed  to  provide  a  means  for  fuel  flow  control 
and  measurement. 

The  models  tested  in  the  small  APL  facility  are  shown  schematically 
in  inverted  position  in  Fig.  26.  Wedge  angles  and  major  dimensions  are  given 
in  the  code  chart.  Side  plates,  to  prevent  transverse  spillover,  were  avail¬ 
able  for  all  models  with  width  (d)  less  than  the  tunnel  width  of  5.92  inches. 
Fuel  was  injected  from  a  number  of  0.031-in. -diam  fuel  ports  on  1/2-in.  span- 
wise  spacing  supplied  by  one  of  several  available  fuel  manifolds  located 

at  or  slightly  upstream  of  the  knee.  Generally,  the  models  were  tested  at 

c  o 

angles  of  attack  from  -4  to  +4  and  with  fuel  flow  rates  from  the  lean 

ignition  limit  to  a  maximum  rate  above  which  thv  propulsion  tunnel  would 

not  operate  due  to  excessive  heat  release. 
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ERJ-W-10-30 


CODE  DESIGNATION 

*1 

S2 

d  (in.) 

r  (in.) 

C  (In.) 

MBS! 

ERJ-  W-10-10 

10° 

10° 

S.92 

1.03 

12.10 

ERJ -W- 10  -20 

10° 

20° 

4.00 

1.01 

8.51 

ERJ -W- 10- 30 

10° 

30° 

4.00 

1.01 

7.44 

4.19 

ERJ -W-(  10-  22  -  38  )  - 10 

I0°,22°,38° 

10° 

4.00 

1.01 

9.03 

Fig.  26  SKETCHES  AND  DIMENSIONS  OF  EXTERNAL  BURNING  MODELS 
TESTED  IN  THE  6  !N.  BY  7  IN.  TUNNEL  AT  APL 
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Typical  axial  and  transverse  static  pressure  distributions  generated 
with  triethyl  aluminum  (TEA)  fuel  on  ERJ-W-10-30  at  +4°  angle-of-attack  are 
shown  in  Fig.  27.  The  pressure  rise  due  to  combustion  on  the  expansion 
wedge  increased  with  increasing  fuel  rates,  w^,  but  the  pressure  rise  per 
lb/sec  of  fuel  dp/dw^)  decreased.  In  this  test  the  lean  ignition  point  was 
near  -  0.009  lb/sec,  but  in  other  tests  with  other  models  at  different 

angles  of  attack,  lean-limit  fuel  rates  as  low  as  0.0036  lb/sec  were  obtained. 

Air  passing  over  the  model  is  compressed  on  the  forward  surface  from  a 

free  stream  static-pressure/total-pressure  ratio  of  0.0018  to  a  ratio  of  0,008 

by  turning  through  the  oblique  bow  shock.  In  the  absence  of  heat  release, 
the  f 3 ow  r<. -expands,  at  the  knee,  to  a  pressure  level  near  that  for  a  simple 
Prandtl-Meyer  turn.  When  the  lean-limit  fuel  flow  is  established,  the  volume 
source  due  to  thermal  expansion  in  the  heat  release  zone  essentially  cancels 
the  Prandtl-Meyer  expansion  in  the  adjacent  (external)  air  streamtube,  so  that 
the  pressure  on  the  expansion  surface  is  nearly  the  same  as  that  on  the  com¬ 
pression  surface.  Increasing  the  fuel  rate  increases  the  strength  of  the 
oblique  shock  at  the  injection  station  and  produces  higher  forces  on  the  aft 
wedge.  Thus,  net  thrust  can  be  produced.  However,  as  previously  noted,  dp/dw^ 

is  negative  at  the  higher  w^  so  that  specific  impulses  for  thrust  and  lift 
decrease  as  w^  is  increased.  This  effect  can  be  seen  in  a  summary  plot  of 

the  data  from  ERJ-W-10-30  shown  in  Fig.  28.  The  data  from  Fig.  27  are  plotted 
as  curve  F,  and  since  the  slope  for  either  AC7  or  AC  versus  fuel  flow  rate 
is  less  than  that  of  the  corresponding  theoretical  curve  for  O'  =  +  4°, 
efficiency  decreases  as  fuel  rate  increases.  The  only  data  with  a  different 
trend,  curves  B  and  C  were  with  fuel  heated  above  the  nominal  level  of  200°  F. 
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The  theoretical  curves  of  Fig.  28  are  based  on  the  constant  pressure 
analysis  of  Section  2.3  (Fig.  15)  using  C  1^  fuel  at  ER  =  1.0  rather  than 
for  the  actual  fuels  used  in  the  tests.  The  effect  of  fuel  change  would  be 
small  because  thh  fuel  flow  required  for  a  given  pressure  rise,  hence,  lift 
and  thrust,  is  primarily  dependent  on  the  heating  value  of  the  fuel.  TEA 
has  a  heating  value  of  18,360  Btu/lb  and  TIBA  (tri-isobutyl  aluminum)has 
18,420  Btu/lb  as  compared  to  a  nominal  value  of  18,630  Btu/lb  for  C  H„  from 
Ref.  31  .  Theoretical  curves  for  ER  lower  than  1.0  would  be  somewhat  better, 
(see  Fig.  16(D))  but  since  no  instream  gas  sampling  or  measurement  of  the 
capture  height  was  attempted,  the  effective  experimental  equivalence  ratio 
cannot  be  determined.  Moreover,  the  theoretical  curves  do  not  consider  non¬ 
equilibrium  kinetic  effects  and/or  mixing  effects,  which  for  small  models 
probably  are  very  influential  in  the  trends  of  performance  with  ER. 

In  general,  the  results  approach  theoretical  lift  gains  more  closely 
than  theoretical  thrust  gains  due  to  the  fact  that  the  fuel  is  injected  just 
ahead  of  the  model  knee.  The  small  force  component  due  to  injection  and 
combustion  on  the  forward  wedge  adds  to  lift  but  subtracts  from  thrust,  whereas 
in  the  theoretical  model  all  pressure  rise  effects  are  felt  only  on  the 
rearward  facing  surface.  The  theoretical  values  are  approached  only  at 
very  low  fuel  rates  or  with  externally  preheated  fuel.  TEA  gives  better 
results  than  TIBA  for  a  given  fuel  rate  at  these  conditions  (compare  A  and 
D). 

Pressures  on  the  aft  plate  extension  were  generally  about  the  same 
as  on  the  rearward  facing  wedge,  indicating  that  heat  release  was  still 
occurring  downstream  of  the  wedge.  This  pressure  increase  would  increase 
lift  at  all  angles  of  attack  but  would  subtract  from  thrust  for  positive  O' 
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and  add  at  negative  O'.  Residence  time  based  on  the  air  velocity  would  be 
about  25  p  sec  in  the  aft  wedge  zone  and  an  additional  60  p  sec  in  the 
extension  plate  zone.  Results  from  the  other  models  tested  were  similar, 
with  most  efficiencies  between  25  and  60  percent  and  a  few  very  lean  operating 
points  having  about  the  theoretical  efficiency. 

The  larger  model  (Fig.  29)  used  in  the  free- jet  tests  at  OALwas 
geometrically  similar  to  ERJ-W-10-20  tested  at  APL  but  was  three  times 
larger  in  linear  dimensions.  Each  fuel  manifold  fed  a  row  of  twenty-one 
0. 047-in. -diam  injection  holes  on  0.5-in.  spacing.  The  two  manifolds  used 
on  these  tests  were  located  1.60  in.  and  0.25  in.  upstream  of  the  knee  and 
were  designated  F2  and  F3,  respectively.  Figure  29(B)  is  a  photo  of  the 
basic  model  installed  in  the  32-in. -diam,  Mach  5  free- jet  nozzle.  The 
removable  fuel  injection  pylons  shown  in  Fig.  29(C)  have  30°  total-included- 
angle  front  wedges  and  40°-total-included-angle  rear  wedges.  Six  equally 
spaced,  0.047-in. -diam  fuel  ports  (24  total)  were  located  on  each  obtuse 
angle  knee  of  a  pylon.  Each  pylon  provides  its  own  bow  shocks,  which  lie 
transversely  within  the  field  of  the  main  body's  bow  shock.  These  additional 
shocks  double  the  static  pressure  and  raise  the  static  temperature  by  about 
507»,  i.e.,  from  400°R  to  600°R.  Aluminum  oxide  deposit  patterns  can  be 
seen  on  the  surface  in  Fig.  29(C).  The  presence  of  oxides  upstream  of  the 
injector  pylons  indicates  that  a  boundary  layer  separation  existed  during 
combustion. 

The  maximum  air  total  temperature  that  could  be  provided  by  the 
facility  was  1500°F,  (  ~  700°R  below  that  required  for  Mach  5.0  flight 
simulation),  hence  the  free  stream  static  temperature  was  ~  250°R,  or 
about  140°R  low,  and  some  reservation  is  needed  in  judging  scaling  effects 
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A)  SCHEMATIC  MODEL  WITH  AFT  EXTENSION  PLATE  BUT  WITHOUT  SIDE  PLATES 


B)  BASIC  MODEL  INSTALLED  AT  EXIT  OF  MACH  5 
FREE  JET  NOZZLE.  FUEL  LINES  COME  OFF  NEAR 
SIDE,  INSTRUMENTATION  LINES  OFF  FAR  SIDE. 


C)  TOP  REAR  VIEW  OF  MODEL 
WITH  FUEL  INJECTION  PYLONS, 
SIDE  PLATES,  AND  AFT  PLATE 
INSTALLED. 


Fig.  29  ERJ  MODEL  TESTED  AT  ORDNANCE  AEROPHYSICS  LABORATORY 
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between  the  APL  and  OAL  tests-. 

Axial  pressure  distributions  for  tha  model  with  side  plates  and 
aft  plate  extension  for  0°  and  4°  acgle-of-attaek  are  given  in  Fig.  30. 

For  the  former,  (A)  injection  is  from  the  forward  F2  manifold  only,  whereas 
for  the  latter,  (B)  injection  from  F2  and/or  F3  is  used.  In  all  burning 
runs  the  pressure  on  the  expansion  wedge  is  more  or  less  constant.  At  most 
fuel  rates  injection  produced  a  pressure  spike  which  rapidly  dissipated  at,  or 
just  downstream  of,  the  knee.  In  these  tests  (as  in  all  others),  the  pressure 
spike  was  located  farther  aft  with  injection  from  the  F3  manifold.  Some  of 
the  effects  noted  in  tests  with  this  and  other  configurations  and  test 
conditions  are: 

1.  For  the  same  ratio  of  fuel  flow  to  air  total  pressure,  results  were 
better  at  a  total  pressure  of  100  psia  than  at  160,  200  or  70  psia.  If  jet 
penetration  depends  primarily  on  the  relative  fuel-air  momentum,  then  the 
penetration  of  a  incompressible  liquid  into  a  compressible  gas  will  increase 
with  increasing  air  pressure  for  constant  fuel  port  area  and  ratio  of  fuel 
flow  to  air  pressure;  e.g.,  doubling  both  the  fuel  rate  and  the  air  pressure 
quadruples  the  fuel  momentum  but  only  doubles  the  air  mome”tum.  Since  the 
penetration  (or  effective  Y^)  increases,  the  effective  equivalence  ratio 
decreases,  and  it  is  conjectured  that  in  these  tests  the  optimum  ER  condition 
occurred  near  100  psia.  Optimum  ER  again  involves,  not  onl;,  the  optimization 
of  the  equilibrium-one-dimensional  process  per  the  theoretical  analysis  but 
also  includes  the  non-equilibrium  effects  (reaction  rate),  atomization  of 
the  fuel,  mixing  and  the  effects  of  a  hot  air  boundary  layer. 
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o=  0,  M0=  5.0,  T,0=  1550°R,  p,Q=  97  psia,  TEA  FUEL,  F2  MANIFOLD 


A)  AXIAL  PRESSURE  DISTRIBUTION 


a  =  -44,  M0  =  5.0,  T,0  =  i530°R,  Pt()=  98  p*ia,  TEA  FUEL,  F2  AND/OR  F3  MANIFOLD 


B)  AXIAL  PRESSURE  DISTRIBUTION 

Fig.  30  EXPERIMENTAL  PRESSUR  E  PROFI LES  FROM  OAL  TESTS  OF  THE  10°-20° 
WEDGE  MODEL  WITH  SIDE  PLATES  AND  AFT  PLATE. 
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2.  Results  with  the  pylons  were  disappointing  in  that  only  slightly 
higher  pressure  fields  were  produced  as  compared  to  tests  without  pylons. 
Injection  from  the  pylons  only  was  more  effective  than  from  combined  injection 
from  pylons  and  body. 

3.  Combination  schlieren-direct  luminosity  photographs  (with  a  slow 
shutter  on  the  direct  lens)  of  the  basic  model  show  the  presence  of  relatively 
strong  injection  shock  and  weaker  combustion  induced  shocks.  Figure  31(A)  is 
a  typical  photograph  showing  a  luminous  main  flame  zone  and  a  small  zone  of 
burning  in  the  boundary  layer  just  upstream  of  the  fuel  injection  ports. 

In  most  cases,  the  boundary  layer  flame  traveled  3  to  4  inches  upstream  of 
the  injection  point  and  caused  a  weak  oblique  shock  in  the  external  flow. 
Figure  31(B)  shows  a  case  with  poorer  combustion  than  in  Fig.  31(A).  There 
appears  to  be  a  region  of  relatively  low  heat  release  (less  luminosity)  in 
the  first  portion  of  the  flame  zone  and  a  weaker  flame  shock.  Inspection 
of  the  corresponding  pressure  plots  shows  that,  in  this  case,  the  pressure 
rise  due  to  injection  was  followed  by  expansion  around  the  knee  before  the 
secondary  flame  shock  was  reached,  at  which  point  the  pressure  increased 
about  25%  and  remained  constant  to  the  trailing  edge. 

4.  In  general,  the  large-scale  tests  (at  OAL)  showed  lesser  thrust 
and  lift  gains  per  pound  of  injected  fuel  (50-75%)  than  the  smaller-scale 
tests  (at  API.).  Besides  the  aforementioned  lower  temperature  of  the  OAL 
tests,  it  was  suggested  that  dissimilar  (cooler  and  relatively  thinner) 
boundary  layer  and  relative  fuel  jet  penetration  may  have  contributed  to 
the  scaling  effects. 
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B)  INJECTION  FROM  F3,  w,  =  0.20  lb  sec,  p,0  =  195  ps.o,  ct  =  0° 


Fig.  31  SCHLIEREN  PHOTOGRAPHS  FROM  OAL  TESTS  OF  BASIC  WEDGE  MODEL 
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The  flat-plate  model  tested  in  the  third  phase  of  the  program 
is  shown  in  Fig.  32(A).  It  is  a  water-cooied  plate  6  inches  wide,  ]5  inches 
long  and  1.25  inches  thick.  In  the  tests  reported  herein,  fuel  was  injected 
from  the  most  forward  fuel  manifold  which  contained  twelve  0.031  in.  diam. 
equally  spaced  fuel  ports.  The  model  surface  contained  numerous  static 
pressure  taps,  and  a  traversing  pitot  tube  was  used  to  make  surveys  of  the 
flow  field  in  the  combustion  zone.  In  the  tests  described,  the  Mach  number 
in  the  flow  field  ahead  of  the  fuel  injection  shock  was  5.04.  The  static 
pressure  was  0.78  psia  and  the  static  temperature  was  345-365°R,  25-45°R 
lov  for  true  simulation  at  66,000  ft. 

Figure  32(B)  shows  static  pressure  profiles  taken  along  the  model 
centerline  for  "cold"  flow  (no  injection)  and  for  toluene  and  triethyl 
aluminum  (TEA)  fuel  injection  from  the  forward  fuel  manifold.  The  rather 
abrupt  pressure  rise  2.5  in.  from  the  leading  edge  followed  by  a  decay  back 
to  about  the  free-stream  level  in  the  toluene  curve  is  as  expected  for  non¬ 
reacting  liquid  injection.  Schlieren  observations  of  the  flow  field  confirmed 
this  results. 

With  TEA  injection  the  pressure  first  rises,  then  decays  slightly 
due  to  the  injection  shock  followed  by  weak  expansion  waves,  and  then  rises 
to  a  near  constant  value  as  the  heat  release  starts  to  take  place.  The 
continuing  relatively  constant  pressure  plateau  to  the  trailing  wedge 
indicates  that  heat  release  is  continuing  and  has  not  been  completed  in 
the  12  inches. 
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REPLACEABLE  VARIABLE 

LEADING  EDGE  FUEL  INJECTORS 


A)  SCHEMATIC  ILLUSTRATION  OF  MODEL 


DISTANCE  FROM  LEADING  EDGE  (inches) 


B)  LONGITUDINAL  PRESSURE  PROFILES 


Fig.  32  FLAT  PLATE  COMBUSTOR  MODEL  TESTED  AT  MACH  5.04  (Ref.  34) 
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In  order  to  appraise  the  latter  situation  a  study  of  droplet  evapora¬ 
tion  and  diffusion  into  a  second  medium  was  undertaken.  The  technique  pro¬ 
posed  by  Zwick,  Grubman  and  Hardy,'-3  which  considers  spherical  particles 
with  known  drag  and  heat  transfer  coefficients  was  used  to  determine  the 
droplet  diameter  as  a  function  of  time.  A  necessary  input  for  this  calcu¬ 
lation  is  the  initial  droplet  diameter.  The  empirical  relationship  of 
Ingebo  and  Foster-3  ^ 

d0  - 3-9  [vi/Wg3)  °'25  d°' 5  W9) 

was  used  to  determine  the  volume  mean  initial  drop  diameter.  For  the  TEA 
flow  rate  of  9.047  lb/sec  at  230°F  the  calculated  dQ  is  .290  x  10"4ft,  or 
about  9  microns.  Diffusion  coefficients  of  TEA  into  air  were  determined  using 
the  rigid  sphere  model  of  Hirschfelder ,  at  al5 7  based  on  estimated  molecular 
diameters  for  TEA. 

To  calculate  the  droplet  trajectory  the  analysis  by  Zwick,  et  al-3-3 
had  to  be  modified  to  account  for  the  observed  combustion  phenomena.  In 
their  analysis  it  was  assumed  that  the  droplet  is  accelerated  by  an  air  stream 
parallel  to  the  plate  and  therefore  the  fuel  penetration  is  dependent  only 
on  its  initial  momentum  and  the  ratio  of  injection  and  free  stream  velocities 
and  the  physical  properties  of  the  fluids.  The  penetration  resulting  from 
this  type  of  calculation  would  be  but  a  small  fraction  of  an  inch  at  the 
trailing  edge.  This  is  contrary  to  the  observed  wedge-shaped  luminous  region 
extending  from  the  surface  at  two  inches  aft  of  the  injection  station  to 
about  1.8  inches  from  the  plate  at  the  trailing  edge. 


-85- 


TH«  JOHNS  HO^KI HS  UMVt«S»TV 

APPLIED  PHYSICS  LABORATORY 

Sitvf*  &+niHc  ma«h4m? 


i 


i 


Greater  fuel  penetration  can  be  postulated  in  the  following  way. 

The  fuel  .ene-tr  tea  through  the  separated  region  with  no  effective  down¬ 
stream  disj.lacc-.iant  (in  fact,  some  fuel  can  recirculte  upstream  and 
burn  in  the  ce  ".rated  boundary  layer,  as  in  the  OM  toots),  «L  which  point 
It  encounters  r.n  uirsfcr*_-ani  moving  not  parallel  to  the  plate,  but  away  fr  •  , 
the  surface  at  an  inclination  det( rmined  by  the  injection-separation  shock 
strength.  The  air  stream  accelerates  the  droplet  both  downstream  and  away 
from  the  .  late.  Subsequent  turning  cf  the  air  stn.  due  to  the  expansion 
and  re  com.  rear, ion  by  the  flame  shock  similarly  affects  the  droplet,  Gnce 
the  flume  shock  is  cro:sed,  the  direction  of  the  main  stream  4 ..  dependent  on 
the  distance  fr'.m  the  plate  according  to  the  postul  '.ad  constant  pressure 
heat  addition  model.  The  stream  moves  parallel  to  the  plate  near  the  sur¬ 
face  but  in  the  direction  of  the  main  stream  doflecli  .n  along  the  hot-c.ld 
interface. 

Three  calculated  droplet  trajectories  are  shown  in  Fig.  33(d)*  The 
lowest  curve  represents  a  clculation  o.ccordin0  to  the  technique  by 

'"r 

Zwiclc,  ct  al.^''  Ti.o  t'.'c  upper  curves  were  determined  ..ccording  to  the 
above-described  Modification  of  that  ti.eory.  In  the  U\.joctorJ  will,  the 
greatest  j.cnetr  •  tisn  the  droplet  is  assumed  to  be  near  the  hot-cold  inter¬ 
face.  In  the  othei  case,  It  is  assumed  to  be  near  the  center  of  the  heat 
addition  sons .  These  curves  accordingly  diverge  ;:.,'3  i-  •  lownt.tr- am  of  the 
injection  ^olnt  ’..’here  they  enter  the  flame  sons.  For  reference,  both  the 
observable  lu...inou~  zone  and  the  position  of  the  bound 'ry  sf  heat  rele.se 
as  deduced  .  i.tct  .return  measur..i.ientj  aie  shown. 
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A.  OBSERVED  LUMINOUS  BOUNDARY 

B.  HEAT  ZONE  BOUNDARY  FROM  PITOT  MEASUREMENT 

C.  TRAJECTORY  BASED  ON  INTERFACE  DEFLECTION 

D.  TRAJECTORY  BASED  ON  CENTERLINE  DEFLECTION 

E.  TRAJECTORY  BASED  ON  NO  NORMAL  ACCELERATION 


DISTANCE  FROM  INJECTOR  (inches) 


A)  DROPLET  TRAJECTORIES 


MACH  NUMBER 


B)  SCHEMATIC  ILLUSTRATION  OF  FLOW  FIELD 


Fig.  33  SCHEMATIC  ILLUSTRATION  OF  FLAT  PLATE  COMBUSTOR  FLOW  FIELD 
AND  CALCULATED  DROPLET  TRAJECTORIES. 
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In  order  to  analyze  the  pitot  pressure  measurements  it  is  necessary 
to  postulate  a  more  applicable  flow  model  than  that  used  in  the  theoretical 
study  (Fig.  18),  since  combustion  is  displaced  a  slight  distance  downstream 
and  is  only  partially  completed.  A  viscous  layer  is  also  present.  Figure 
33(B)  shows  the  concept  of  the  flow  picture.  The  shock  wave  patterns  and 
flame  zone  picture  are  consistent  with  schlieren  and  luminosity  photographs 
of  the  flow  field. 

In  this  model  the  incoming  air  is  divided  into  three  major  streamtubes: 
boundary  layer  air,  inviscid  burned  air,  and  air  that  is  turned  but  not  burned. 
The  typical  Mach  number  profile  indicates  that  the  properties  are  not  uni¬ 
form  in  the  major  streamtubes  as  is  assumed  to  be  the  case  in  one-dimensional 
sttf8ies.  Thus,  for  detailed  analysis  the  streamtubes  must  be  further  sub¬ 
divided  or  suitable  integrations  of  properties  must  be  made  to  reduce  the 
problem  to  a  one-dimensional  analysis.  The  major  streamtube  boundaries  are 
determined  in  the  following  manner:  The  height  of  the  boundary  layer  stream- 
tube  at  the  injection  shock  h^  is  determined  by  a  cold-flow  pitot-pressure 
traverse  and  checked  with  a  calculated  boundary  layer  thickness.  This 
calculation  also  provides  an  estimate  for  the  total  mass  flow  in  the  streamtube 
before  fuel  injection.  After  fuel  injection  it  is  assumed  that  no  further 
flow  into  the  boundary  layer  occurs  and  that  the  increase  in  h^  above  h^ 
is  due  to  heat  release  and  additional  viscous  losses.  The  downstream  pitot 
traverse  determines  h,  and  the  sum  of  h„  and  h, •  Since  the  flow  in  the 
outside  streamtube  is  inviscid  and  adiabatic,  the  upstream  height  h^  can 
be  specified  by  solving  wedge-flow  relationships  consistent  with  observed 
shock  waves  and  pressure  measurements.  The  difference  between  the  total 
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height  of  the  three  streamtubes  and  h^  +  h^  determines  h Assuming  that 

the  injected  fuel  is  proportionally  divided  on  a  mass  basis  between  the  two 

major  heat  release  streamtubes,  then  all  of  the  necessary  conditions  are 

now  specified  to  analyze  the  combustion  except  the  J  pdA  term  in  the  momentum 

34 

equation  for  each  streamtube.  This  integral  is  the  product  of  the  local 
pressure  and  the  incremental  projected  area  in  the  flow  direction  summed 
over  the  streamtube  boundaries. 

As  a  first  approximation,  the  wall  static  pressure  is  assumed  to  be 
constant  in  planes  normal  to  the  surface.  A  further  refinement  can  be  made 
by  analyzing  schlieren  pictures  of  the  flow  field  and  actually  positioning 
the  shock  and  expansion  waves  and  calculating  the  static  pressure  field. 

If  the  Mach  profile  across  h,.  is  relatively  flat,  then  the  center  streamtube 
can  be  individually  handled  as  a  one-dimensional  flow  process.  In  this 
circumstance  the  pressure-area  integral  term  is  found  by  assuming  the  wall 
static  pressure  on  the  plate  side  of  the  streamtube  and  using  the  wedge 
flow  pressures  on  the  outer  side. 

For  the  data  point  shown,  the  TEA  fuel  flow  rate  was  0.049  lb/sec. 

The  pitot  profile  measurements  showed  h,.  +  h^  =  1.20  in.  The  calculated 
values  for  h^,  and  h^  were  0.05,  0.61,  and  1.93  in.,  respectively.  At 
the  traversing  plane,  which  was  7.3  in.  downstream  of  the  injection  porl , 
the  calculated  average  combustion  efficiency  was  31.5%.  The  average  equiva¬ 
lence  ratio  in  the  combustion  zone  was  almost  1.0. 

It  is  possible  to  relate  the  experimental  result  to  the  more  elementary 
theoretical  model  of  Fig.  18.  For  the  calculated  and  the  measured  pres¬ 
sure  ratio  of  0.0038/0.0015  =  2.53  at  Mach  5.04,  the  theoretical  length  of 
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the  constant  pressure  region,  X^,  is  31.0  in.,  and  the  total  length  is  62. A  in. 
Combustion  would  be  completed  24  in.  from  the  leading  edge  and  the  resulting 
lift  specific  impulse  would  be  5760  sec.  An  integration  of  the  pressure  rise 
in  the  test  r  ported  with  the  measured  fuel  rate  gives  a  normal  force  specific 
impulse  of  1350  sec. 

From  the  theoretical  constant-pressure  heat  addition  model  it  is  also 
possible  to  show  the  variation  of  total  temperature  with  distance  in  the  heat 
release  zone  for  a  given  capture  height,  Y^.  Thus,  for  the  constant-pressure 
process  the  local  velocity  and  Mach  number  are  known  as  a  function  of  distance. 
In  tests  including  the  above  made  at  nearly  the  same  condition  extensive 
pitot  pressure  surveys  were  made  and  the  local  Mach  number  was  deduced.  The 
Mach  numbers  based  on  measurements  made  along  the  centerline  streamtube  of 
the  heat  release  zone  are  shown  in  Fig.  34  and  are  compared  to  the  Mach  number 
variation  with  distance  from  the  theoretical  analysis.  The  distance  for  the 
experimental  results  is  the  distance  along  the  streamline  from  the  point  of 
first  perceptible  heat  release  to  the  pitot  measuring  station.  The  close 
correlation  of  data  with  experiment  suggests  that  the  theoretical  description 
of  the  heat  release,  viz,  a  continuous  rise  in  total  temperature  with  distance, 
is  reasonable  for  these  tests. 

58  o 

Hypersonic  wind  tunnel  tests  were  made  at  Boeing  with  the  8  half¬ 
angle  blunted  cone  model  shown  in  Fig.  35 .  Fuel  was  injected  from  a  single 
orifice  located  at  the  one-quarter  body  length  point.  The  fuel  orifice  diameter 
was  varied  from  0.003-in.  to  0.024-in.  and  the  orientation  was  either  normal 
to  the  surface  or  inclined  45°  upstream.  Although  the  model  contained  pres¬ 
sure  instrumentation  and  was  mounted  in  a  two  component  force  balance  (normal 
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force  and  pitching  moment),  both  measurements  were  not  made  simultaneously. 

In  a  typical  run  the  tunnel  flow  was  stabilized,  the  sting-mounted  model 
was  inserted  into  the  air  stream,  fuel  flow  was  established,  and  the  model 
was  rotated  from  -30°  to  440°  angle-of-att  Either  pressure  or  force 

data  were  taken,  and  schlieren  and  direct  luminosity  photographs  were  made. 

The  tests  were  made  in  a  12-in. -diameter  open- jet  tunnel  which 
was  operated  at  a  maximum  stagnation  pressure  of  1000  psia  and  maximum 
stagnation  temperature  of  1460°R.  Following  some  preliminary  cold  flow 
and  nonreacting  injection  runs  a  single  test  at  Mq  =  6.1  was  made  using 
pentaborane  injected  fror.:  a  0.012-in. -diameter  hole  oriented  normal  to  the 
surface.  The  free  stream  static  conditions  for  this  tests  were  0.4  psia 
and  140°R,  or  250°R  low  for  true  simulation.  With  a  fuel  flow  rate  of 
0.0048  lb/sec  and  with  the  model  at  zero  angle-of-attack,  no  ignition  of  the 
fuel  was  observed  even  with  a  spark  ignition  source  located  at  the  base 
of  the  model.  Further  testing  at  this  Mach  number  was  abondoned  and  instead 
a  Mach  5.0  nozzle  was  used  in  the  remainder  of  the  tests.  All  combustion 
tests  at  this  Mach  number  were  made  with  free  stream  static  conditions  of 
0.75  psir  and  250°R,  or  140°R  low  for  simulation.  In  the  first  group  of 
successful  tests  at  this  condition,  pentaborane  was  injected  from  0.009- 
and  0.012-in. -diameter  orifices  at  45°  upstream  and  normal  to  the  surface. 
Direct  luminosity  photographs  revealed  that  no  combustion  occured  at  angles 
of  attack  from  -30°  to  30°.  (Again,  o'  is  positive  for  windward  side  burning.) 
However,  a  combustion  flame  was  established  on  the  aft  portion  of  the  model 
at  o'  =  +40^  with  all  fuel  injector  configurations.  Figure  36  shows  the 
static  pressure  traces  for  the  different  injector  conditions  tested  at 
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AXIAL  LOCATION  (inch.*) 


Fig.  36  AXIAL  STATIC  PRESSURE  DISTRIBUTION  FOR 
CONE  TEST  WITH  PENTABORANE  AT  a  =  40° 
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O'  =  440°.  In  all  cases  a  region  of  negative  pressure  coefficient  occurs 
downstream  of  the  injection  location.  Presumably  this  is  the  overexpansion 
region  following  a  pressure  spike  at  the  injector,  which  was  not  detected 
in  the  relatively  widely  spaced  pressure  taps.  In  support  of  this  conjecture 
are  both  the  schlieren  pictures  which  show  the  presence  of  a  strong  oblique 
wave  at  the  injector  and  the  higher  pressure  reading  on  tap  4  in  run  16. 

No  explanation  of  the  high  readings  on  tap  2  in  runs  16  and  20  is  given. 

The  pressure  begins  to  rise  approximately  one  inch  downstream  from  the 
injector  due  to  the  observable  heat  release  and  remains  positive  for 
approximately  4  inches  farther  downstream  and  then  drops  due  to  severe 
wrap-around  at  these  high  angles  of  attack.  Considering  the  discrepancy  between 
the  data  from  runs  14  and  20  at  essentially  identical  test  conditions,  no 
conclusions  can  be  made  regarding  the  effects  of  fuel  orifice  diameter  and 
angular  orientation  on  performance.  The  base  ignitor  was  used  in  all  of 
these  tests  and  when  damage  was  sustained  by  it  in  run  20,  a  switch  to  more 
reactive  aluminum  borohydride  was  made.  The  ignition  source  was 

not  used  in  the  remainder  of  the  tests. 

Using  45°  upstream  injection  from  a  0.012-in. -diameter  hole,  two 
runs  were  made  with  Aj^BH^i  the  first  a  pressure-measuring  run  with  a  from 
-30°  to  440°  and  the  second  a  force-balance  run  with  a  from  -25°  and  4-25°. 

The  fuel  flow  rates  were  0,0066  lb/sec  and  an  estimated  0.0076  lb/sec 
respectively.  A  series  of  direct  luminosity  photographs  are  shown  in 
Fig.  37  for  a  =  -30°  to  440°.  Combustion  initiated  at  the  fuel  orifice 
for  all  model  orientations  and  contined  downstream.  With  windward 
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a=  0° 


a  =  -30° 


Fig.  37  LUMINOSITY  PHOTOGRAPHS  OF  ALUMINUM 

BOROHYDRIDE  COMBUSTION  ON  AN  8°  CONE  (Ref.  58) 
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burning  (positive  O')  severe  wrap-around  is  observable  at  a  =  40°  and  to  a 
lesser  degree  at  a  =  20°.  Pressure  traces  for  a  =  0°,  +20  and  440°  are 
shown  in  Fig.  38;  plugging  of  the  pressure  taps  with  solid  deposits  pre¬ 
vented  pressure  measurements  at  negative  a-  Schlieren  photographs  indicate 
a  significantly  stronger  disturbance  at  the  injector  than  is  implied  by 
linear  interpolation  from  tap  4  to  tap  5.  The  general  trend  of  the  pressure 
traces  is  a  pressure  decay  following  the  initial  injection  compression  and 
then  a  pressure  rise  due  to  heat  release  which  is  significantly  different 

for  the  three  data  runs.  The  writer  made  an  estimate  of  C„  based  on  the 

N 

pressure  readings  as  follows.  Assuming  that  the  measured  pressures 
applied  to  the  longitudinal  projected  area  of  the  model  (since  the  luminous 
flame  appeared  to  have  about  the  same  width  as  the  model  at  any  station), 
was  computed  by  integrating  the  pressure  profile  over  this  projected  area. 

At  a'  =  +20°  and  +40°  a  negative  force  coefficient  was  obtained  and  at  a  =  0° 
a  positive  value  of  CN  =  0.13  was  computed.  The  force  balance  runs  made 
at  similar  test  conditions  were  interpreted  (transient  heating  of  the  balance 
caused  significant  zero  shifts)  to  give  a  =  0.035  at  zero  angle-of-attack 

and  CN  =  0.03  at  a  ~  -25°.  The  normal  force  changed  sign  at  approximately 
+3°  angle-of-attack  and  dropped  0.06  below  the  r.on-burning  value  at  o'  =  +25°. 
Specific  impulses  for  the  two  tests  at  a  =  0  would  be  .13  x  68. 8/. 0066  = 

1355  sec  based  on  the  integrated  pressure  run  and  0.035  x  68. 8/. 0076  =  316  sec 
from  the  force  balance  run. 

There  was  virtually  no  change  in  pitching  moment  data  at  negative 
and  zero  a  which  points  out  the  difficulty  in  obtaining  attitude  control 
with  forebody  injection  when  the  effects  of  heat  release  are  being  felt  both 


-97- 


i>t  »><i  *o***rt  t*^r*t*yrr 

Atnno  physics  ia»o«ato«y 

Volt  1M*C 


0 


AXIAL  LOCATION  (inch) 

3  4  5  6  7  8 


Fig.  38  AXIAL  PRESSURE  DISTRIBUTION  FOR  CONE  TESTS  WITH 
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fore  and  aft  of  the  model  c.g.  The  change  in  moment  coefficient  (Cw)  was 

sizeable  at  positive  a,  increasing  from  about  0  at  a  -  0°  to  ~  0.025  at 

o  o 

a  =  16  and  then  rapidly  decaying  to  0  at  O'  =  25  .  This  strange  behavior 

is  due  to  the  large  change  in  the  character  of  the  pressure  traces  at  positive 
O'  and  would  lead  to  control  design  problems  in  a  flight  vehicle  because  the 
nature  of  the  pressure  variations  due  to  combustion  would  probably  be 
strongly  influenced  by  flight  environment,  i.e.,  at  low  altitude,  higher 
pressure  would  accelerate  the  heat  release  rates. 

A  mixture  of  25%  aluminum  borohydride  and  75%  pentaborane  was 
used  in  two  runs  with  the  same  injector  used  in  the  100%.  AX(BH^),  runs.  At 
a  =  0°  with  both  a  very  low  flow  rate  (not  determined)  and  a  flow  rate  of 
0.0052  lb/sec  the  only  combustion  observed  was  in  the  base  region. 

Measurements  were  made  from  the  schlieren  photographs  to  determine 
the  penetration  depth  of  the  fuel  and  comparisons  were  made  with  theoretical 
estimates  from  a  fuel  trajectory  analysis.  Correlation  of  data  with  theory 
for  the  penetration,  defined  as  the  point  where  the  fuel  trajectory  is  inclined 
9°  to  the  airstream  direction,  was  good  for  both  the  non-reacting  injectants 
and  for  the  pentaborane  which  had  not  ignited  at  this  point  on  the  trajectory. 
With  Ai(BH^  the  theoretical  penetration  height  was  similar  to  that  for  penta¬ 
borane,  but  the  experimentally  determined  values  were  approximately  twice  as 
large.  Thus,  combustion  appears  to  significantly  increase  the  fuel  penetra¬ 
tion,  which  is  consistent  with  the  measurements  and  conclusions  of  Ref.  36. 
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4.  SUMMARY  AND  CONCLUSIONS 


Although  conclusions  have  been  drawn  from  the  results  of  particular 
theoretical  and  experimental  studies  in  various  sections  of  this  review, 

it  is  nonetheless  useful  to  consider  the  material  in  its  entirety  with  re¬ 
spect  to  the  state-of-the-art  of  external  burning.  Of  course,  most  of  the 
experimental  work  has  been  exploratory  in  nature  and  more  often  than  not 
was  obtained  at  testing  conditions  less  favorable  for  combustion  than  would 
be  expected  in  flight.  Nevertheless,  a  careful  scrutiny  of  the  available 
data  should  provide  the  systems  designer  with  reasonable  estimates  of  the 
probable  performance  realizable  from  external  burning  and  should  serve  as 
a  guide  to  the  experimentalist  in  planning  additional  tests  which  can  sub¬ 
stantiate  or  negate  the  conclusions  drawn  herein. 

The  experimental  data  are  summarized  in  Table  VII.  The  testing  con¬ 
ditions  and  fuel  flow  rates  are  given  along  with  the  principal  performance 
indices;  normal  and  axial  specific  impulse,  1^  and  1^,  the  maximum 

pressure  coefficient,  C  ,  and  the  normal  force  or  average  pressure  co- 

lmax 

efficient  C^.  The  normal  force  specific  impulse  values  are  plotted  in 

Fig.  39  as  a  function  of  free  stream  Mach  number.  Note  the  break  in  the 

Mach  number  scale  at  Mq  =  5.0  to  accommodate  the  preponderance  of  Mach  5 

data.  For  comparison  two  theoretical  curves  are  shown  for  TEA  combustion, 

one  based  on  linear  theory  and  the  other  on  the  constant  pressure  combustion 

34 

analysis  including  expansion.  The  latter  is  for  the  case  of  a  turning 
angle  equal  to  the  turbulent  separation  value.  The  linear  theory  expression 
can  be  developed  from: 

XfN  =  °p  A  VWf  (50) 
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using  Eqs-  (7)  and  (10) 

«  -  <W>  We  Tc  x£  ; 


C  -  2  6/(X„2  -  1)“ 


A  =  area  under  flame  =  X.d 

i 


’o  '  *  Vo 


0  =  h  w£/w 
f  a 


substituting  Eqs.  (5i)  -  (52)  into  (50), 

"  h  V'  cp  To  «o2  -»4 


=  «  Vo  Y  15 


with  a,  =  V  RT_,  and  S  /  c  -  (Y  ~  1)/Y, 

p 

IfB  -  778  »  <Y  -  I>  W0/a0  <K(,  -  1)* 


(51) 

(52) 

(53) 

(54) 

(55) 


To  obtain  theoretical  values  for  the  other  fuels  used,  the  TEA  values 
can  be  scaled  by  the  ratio  of  heating  value  of  the  particular  fuel  to  that 
of  TEA,  i.e.,  29,360/18,360  =  1.6  for  B_H9,  1.27  for  Ai(BH4)3>  2.81  for  H2 
and  1.0  for  TIEA. 

Nearly  all  of  the  data  fall  considerably  below  the  theoretical 
curves  with  the  notable  exception  of  the  early  NASA  tests'*^  using  ths  very 
reactive  fuel,  Aj£(BH,  )_,  on  relatively  long  models.  Although  for  many  of 

•4  3 

the  tests  it  may  be  argued  that  non-optimum  fuel  infection  systems  produced 
poor  fuel-air  distributions  and  performance  suffered  sccordingly,  it  does 
not  seem  reasonable  to  attribute  all  of  the  deficiency  to  too  rich  or  too 
lean  mixtures,  because  so  many  different  injection  schemes  and  fuel  flow 
rates  have  been  tried.  Rather  it  is  just  as  likely  that  the  limiting 
factor  is  insufficiently  rapid  kinetic,  rates;  i.e.,  for  the  local  condi¬ 
tions  there  was  not  sufficient  time  (distance)  to  complete  the  heat 
release.  This  argument  is  supported  by  the  NASA  results'* \  which  show  a 


-103- 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIFD  PHYSICS  LABORATORY 

SILV2A  SPRING  Makylano 


significant  increase  in  1^  from  run  D  to  run  E  with  an  increase  in 
length  from  21.3  to  46.5-in.,  and  by  the  instream  measurements  made  at 
APL  (Fig.  34),  which  showed  that  the  reaction  was  proceeding  in  the 
streamwise  direction  with  a  nearly  linear  temperature  rise  with  length. 
Similar  effects  were  noted  in  data  from  J  and  K.  On  the  other  hand,  if 
these  were  strictly  fuel-air  distribution  difficulties  one  would  expect 
a  rapid  heat  release  until  all  of  the  local  oxygen  is  consumed,  followed 
by  a  very  slow  release  as  additional  air  mixes  with  the  excess  fuel  and 
products.  Granting  that  this  may  be  a  reaction-rate-limited  situation 
the  question  is:  How  much  length  will  be  needed  to  complete  the  reaction 
in  flight?  Obviously  the  answer  depends  on  the  reactivity  of  the  fuel 
and  on  the  flight  environment.  Aluminum  borohydride  appears  to  be  the 
fastest  reacting  fuel  tested,  followed  by  pentaborane  and  the  aluminum 
alkyls  TEA  and  TIBA.  Higher  pressure  generally  increases  reaction  rates, 
therefore  lower  altitude  should  require  less  length.  Compression  prior 
to  injection  due  to  the  body  shape  and/or  to  burning  on  the  windward  side 
at  angle-of-attack  wi ' 1  help  because  the  local  pressure  and  temperature 
will  be  higher  and  the  velocity  somewhat  lower.  Lower  flight  velocity 
increases  the  residence  time  for  a  given  length  but  also  has  associated 
with  it  lower  maximum  temperatures  in  the  boundary  layer  which  may  be 
important.  Finally,  the  wall  temperature  and  the  fuel  temperature,  which 
could  be  elevated  if  the  fuel  is  used  as  a  regenerative  coolant,  may  also 
be  important-  All  of  these  fa  tors  will  affect  the  answer  to  the  question 
posed;  however;  it  appears  that  the  length  required  to  complete  the  reaction 
will  be  in  the  order  of  feet  rather  than  inches,  judging  from  the  results 
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from  runs  D  and  E  (which  might  be  extrapolated  to  about  6  ft)  and  from 
the  estimate  given  before  of  2  ft  for  reaction  plus  another  3  ft  for  expan¬ 
sion  for  test  M.  If  lengths  of  this  order  are  not  available,  specific 
impulses  are  apt  to  be  lower,  roughly  in  proportion  to  the  ratio  of  the 
available  length  to  the  reaction  length.  Clearly,  there  is  a  need  for 
reaction  rate  data  for  these  fuels  in  order  to  make  better  estimates.  The 
data  possibly  could  be  obtained  in  much  simpler  subsonic  tests  at  simulated 
temperatures  and  pressures. 

For  nearly  all  of  the  data  summarized  in  Table  VII  and  Fig.  39,  the 
length  available  for  combustion  was  quite  short  and/or  the  static  tempera¬ 
ture  was  lower  than  for  simulated  flight,  therefore,  the  impulses  are 
correspondingly  low.  However,  there  are  some  case,  i.e..,.B-E  and  G 
having  roughly  half  the  theoretical  1^,  which  is  much  higher  than  can  be 
obtained  with  non-reactive  thrust  vector  control  systems.  The  perform¬ 
ances  of  a  body  of  revolution, N  ,  was  poorer  than  for  the  two-dimensional 
models,  which  appears  to  be  due  to  excessive  lateral  spillover  and  again 
points  to  the  desirability  of  using  longitudinal  fences.  Thrust  specific 
impulses  were  calculated  for  most  of  the  thrust-generating  configurations 
and  are  listed  in  Table  VII.  Nearly  all  are  low  compered  to  either  theory 
or  competitive  propulsion  systems.  The  highest  value  of  818  sec  for  test 
G  was  obtained  at  the  lowest  fuel  flow  rate,  0.0025  lb/sec.  Conclusions 
bat  d  on  the  theoretical  analysis  regarding  the  performance  potentials  of 
external  burning  systems  in  the  thrust  and  lift  generating  modes,  Figs.  2B 
and  C,  are  adequately  given  in  Section  2.3. 
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The  other  performance  characteristics,  especially  important  in  the 
case  of  the  attitude  controller  (Fig.  2A)  are  the  normal  force  and  maximum 
pressure  coefficients.  For  a  given  normal-force  requirement,  the  control 
surface,  area  varies  inversely  with  the  force  coefficient.  In  most  cases, 
it  is  desirable  to  have  maximum  which,  at  first  glance,  would  seem  to 
be  associated  with  an  external  burner  which  creates  a  very  strong  shock, 
possibly  a  normal  shock.  However,  this  flow  situation  cannot  exist,  because 
at  pressure  ratios  considerably  below  the  normal-shock  value,  the  boundary 
layer  will  separate  ahead  of  the  heat  release  zone  and  the  resulting  pres¬ 
sure  coefficient  due  to  the  volume  source  created  by  the  heat  release  will 
be  more  like  that  of  a  separated  flow.  This  point  is  demonstrated  in 
Fig.  40,  in  which  the  maximum  pressure  and  normal  force  coefficients  for 
the  data  of  Table  VII  are  plotted  against  free  stream  Mach  number.  The  maxi 
mum  pressure  coefficients  occurring  anywhere  in  the  combustion  zone  are 
shown  as  open  symbols  and  the  average  coefficients  in  the  heat  addition 
region  as  closed  symbols.  Most  of  the  data  fall  near  or  below  the  values 
for  turbulent  separation  based  on  Ref.  36  .  This  is  not  to  imply  that 

there  is  a  unique  value  of  C  for  turbulent  separation,  because  there  is 

^max 

some  Reynolds  number  dependence,  but  the  theory  does  represent  a  reasonable 

estimate.  Moreover,  a  lot  of  the  data  were  obtained  at  low  values  of 

Reynolds  numbers  where  the  boundary  layer  should  have  been  laminar,  which 

probably  accounts  for  some  cf  the  lower  C  values  because  laminar  sep- 

59  ^max 

aration  occurs  at  lower  C  .  Even  in  the  few  runs,  B-E  in  which  the 

P 
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pressure  distribution  due  to  injection  and  heat  release  was  characterized 
by  a  pressure  spike  followed  by  a  rapid  decay,  the  maximum  pressure  coeffi¬ 
cients  in  the  spike  region  still  were  of  the  same  order  as  the  turbulent 
separation  value.  Also  shown  for  reference  is  the  variation  for  a 
turning  angle  of  15°,  which  might  represent  a  reasonable  upper  bound  for 
theoretical  system  studies.  Note  that  the  data  in  Fig.  40  were  plotted 
at  the  free  stream  Mach  number  condition  rather  than  at  the  lower  local 
Mach  number  condition  which  would  exist  at  the  injection  point  for  all  but 
the  flat  plate  tests.  Incidentally,  the  normal  shock  value  for  increases 
from  1.250  at  M  =  2  to  1.627  at  M  =  6.5.  These  values  would  be  obtained  if 
the  normal  strong  detonation  (Fig.  6B)  could  be  generated,  but  this  does  not 
seem  possible  when  a  boundary  layer  is  present.  It  appears,  then,  that  the 
designer  of  an  external-burning  system  for  attitude  control  will  have  to 
contend  with  relatively  low  and  CN>  probably  close  to  the  separation  values 
for  the  local  condition  at  injection.  If  higher  force  coefficients  are  required, 
a  combination  of  a  compression  surface  plus  external  burning  will  be  needed. 

Another  consideration  in  the  design  of  external  burning  systems 
which  are  to  provide  iazeral  maneuver  capability  is  the  control  of  the 
moment  coefficient.  With  external  burning  two  choices  seem  plausible.  Either 
the  external  burning  can  be  used  to  provide  the  moment  to  turn  or  trim  the  body 
to  angle-of -attack  to  obtain  the  lateral  force,  or  the  external  burning  can 
provide  the  lateral  force  without  turning  the  body.  In  the  former  case  for  a 
stable  body  it  appears  that  the  fuel  injection  should  be  aft  of  the  body  c.g. 
on  the  lee  side,  rather  than  ahead  of  the  c.g.  on  the  windward  side,  because 
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without  spillover,  toe  heat  release  will  produce  a  positive  pressure  field 
which  will  persist  in  the  streamwise  direction.  Conversely  if  the  external 
burning  is  used  to  trim  an  aerodynamical ly  unstable  body  then  the  injection 
should  be  aft  of  the  body  c.g.  on  the  windward  side.  Thus,  only  for  extremely 
long  bodies  would  the  pressure  field  be  sufficiently  dissipated  with  fore¬ 
body  injection  to  assure  a  positive  moment  coefficient.  Of  course,  one 
might  try  to  exploit  the  wrap-around  characteristics  and  thus  obtain  a 
positive  pressure  coefficient  both  on  the  windward  side  ahead  of  the  c.g. 
and  on  the  reward  side  aft  of  the  c.g.  with  forebody  injection  on  a  stable 
body.  However,  this  scheme  would  seem  to  have  a  very  unpredictable  behavior, 

because  the  pressure  field  would  be  quite  different  at  each  or,  as  was  demon- 

58 

strated  in  the  Boeing  tests. 

For  the  very  rapid  vehicle  maneuvers,  the  non-turning  body  may  be 
attractive.  In  this  case  the  body  remains  at  zero  incidence  to  the  stream. 
Fuel  is  added  on  the  forebody  ahead  of  the  c.g.  on  the  side  opposite  the 
desired  lateral  displacement  in  order  to  produce  a  positive  force  field 
that  is  essentially  balanced,  in  axial  distribution,  about  the  c.g.,  thus 
producing  no  turning  moment.  Fences  of  course  would  be  required.  If  the 
response  time  were  such  that  the  body  could  be  rotated,  then  the  body 
would  be  pitched  as  well,  using  an  auxiliary  control  to  obtain  a  combined 
effect.  The  additional  control  could  possibly  be  an  aft  body  injector  on 
the  opposite  side.  All  of  these  systems  will  require  a  trimming  mechanism, 
which  could  be  a  system  of  small  afterbody  injectors  located  at  varying 
circumferential  positions. 
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A  final  consideration,  not  previoi ^  y  discussed,  is  the  transient 
behavior  of  the  external  burning  system.  For  a  very  rapid  response  system 
using  short  bursts  of  external  burning  for  control,  the  transient  behavior 
of  the  pressure  field  will  have  to  be  known.  To  date  neither  experimental 
nor  theoretical  studies  have  been  made  which  properly  evaluate  the  trans¬ 
ient  behavior. 
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